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ABSTRACT 

An  analytical  model  was  constructed  to  predict  the 
response  of  the  ball  obturator  in  a  ball  obturated,  spin 
stabilized  tubular  projectile  (BOT) .   Torques  acting  on 
the  obturator  were  considered  to  be  due  to  fluid  and  sliding 
friction  forces.   An  experimental  apparatus  was  designed, 
and  built  to  simulate  the  BOT.   Data  obtained  were  the  times 
required  for  the  obturator  to  nutate  through  a  known  angle 
over  a  range  of  projectile  angular  velocities.   Utilizing 
published  coefficient  of  sliding  friction  information  and  a 
linear  approximation  of  the  exact  solution,  the  trend  in  the 
data  to  go  through  a  minimum  value  was  duplicated.   A  least 
squares  fit  of  the  linear  approximation  to  the  data  allowed 
the  forroulation  of  a  new  sliding  friction  coefficient,  y^. 
Using  this  value  of  \i^,    a  plot  of  the  least  squares  fit, 
values  from  the  exact  solution  and  the  experimental  data 
was  constructed.   The  correlation  is  extremely  good  so  that 
confidence  in  the  predictive  capabilities  of  the  model  was 
gained. 
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NOrffiNCLATURE 

(Only  primary  symbols  are  listed.   Intermediate  quantities 

are  defined  in  the  text.) 

A       major  mass  moment  in  inertia,  N-m-s 

C       minor  mass  moment  of  inertia,  N-m-s 

G       applied  torque,  N-m 

h       width  of  gap  between  ball  and  projectile,  m 

M       magnitude  of  applied  torque,  N-m 

r       radius  of  hole  through  projectile,  m 

rg      radius  of  hole  through  ball,  m 

R       radius  of  ball,  m 
s 

S  sin"^  (rs/Rs^ 

9,(}),ijj  Euler  angles  (see  Fig.  7) 

X  (C/A)-l 

y^  fluid  viscosity  (Pa-s) 

Pg  coefficient  of  sliding  friction,  dimensionless 

CO  angular  velocity 

Subscripts 

cr      refers  to  condition  when  Q=Bqj-=^q 

o  initial  value 

p       projectile 

r       refers  to  relative  angular  velocity 

s       sphere  (ball  obturator) 
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x,y,z    refers  to  inertial  coordinate  system 
1,2,3    refers  to  body-fixed  coordinate  system 

Superscripts 

•       denotes  differential  with  respect  to  time 
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I.  INTRODUCTION 

A.   BACKGROUND 

The  tubular  projectile  is  not  a  new  concept.   As  early  as 
1858,  the  usefulness  of  tubular  projectiles  was  recognized  for 
specific  purposes.   In  that  year,  Joseph  Whitworth  (Whitworth 
Threads)  wrote  about  and  included  an  illustration  of  the  tubular 
projectile  in  the  section  on  Rifled  Firearms  of  his  Miscellaneous 
Papers  on  Mechanical  Subjects  [1]^.   The  projectile  pictured 
was  hexagonal  with  a  circular  hole.   Whitworth  noted  its  par- 
ticular "...effectiveness  in  perforating  elastic  materials 
which  prevented  them  from  closing  up."   The  projectile  utilized 
a  wooden  sabot.   Whitworth  also  said  that  the  tubular  projectile 
penetrated  deeper  into  masonry  than  any  he  was  acquainted  with. 

The  Krnka-Hebler  projectile   was  the  next  projectile  refer- 
enced chronologically  [2].   According  to  Charters  and  Thomas, 
this  projectile  was  reported  in  the  Allgemeine  Schweitzerische 
Militarzeitung  as  having  been  so  successful  that  the  Ordnance 
Department  in  the  United  States  carried  out  firing  tests  of 
caliber  .30  tubular  projectiles  in  1894.   A  description  of 
those  test  firings  and  an  updating  of  the  results  was  performed 
by  I.E.  Segal  and  included  in  Ref.  2.   In  the  1894  tests,  a 
vulcanized  fiber  sabot  was  used  to  push  the  projectile  out  the 
barrel. 


Numbers  appearing  in  brackets  refer  to  the  list  of  references. 
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The  purpose  of  the  1894  test  was  to  compare  trajectories 
of  standard  and  tubular  projectiles.   This  was  done  by  com- 
paring the  vertical  drop  on  target  at  a  given  range. 

Segal's  report  in  Ref.  2  indicates  that  the  drag  co- 
efficients computed  from  the  1894  results  agree  closely  with 
those  of  Charters  and  Thomas  [2],  even  though  their  report 
was  published  fifty  years  later  (1944)  and  their  experimental 
apparatus  was  more  sophisticated. 

B.   RECENT  DEVELOPMENTAL  STUDIES 

The  results  obtained  by  independent  researchers  as  to  the 
value  of  the  tubular  projectile  as  an  alternative  to  the  con- 
ventional round  are  contradictory.   It  seems  as  if  there  was 
as  much  research  devoted  to  discrediting  the  tubular  pro- 
jectile as  there  was  to  investigating  its  merits.   Frank  and 
McLaughlin  [3]  have  accumulated  a  great  deal  of  data  from 
various  sources  in  an  attempt  to  "objectively  compare"  the 
merits  of  the  tubular  and  "conventional  shapes".   The  authors 
concluded  that  the  tubular  has  no  particular  advantages  over 
"well  designed"  conventional  projectile  shapes. 

Their  findings  disagree  with  most  of  those  discussed  in 
the  following  sections. 

1.   Range  Tests 

Winchenbach,  Daniel  and  Edgar  [4]  conducted  range 
tests  of  six  configurations  of  tubular  projectiles  and  con- 
cluded that  the  drag  coefficients  were  significantly  lower 
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than  the  standard  High  Explosive  Incendiary  (HEI)  projectile 
of  the  same  caliber.   Only  projectiles  of  the  same  bore  size 
were  compared.   The  experimental  models  were  constructed  from 
standard  20-mm  ammunition  by  boring  various  hole  sizes  to 
obtain  "...  area  ratios  (A^/A^)  of  0.7,  0.8,  0.9  and  1.0." 
Where  A-j-  is  the  throat  area  and  Aj_  is  the  inlet  area  of  the 
projectile.   The  drag  coefficients  of  the  projectiles  with 
area  ratios  of  0.8,  0.9  and  1.0  were  less  than  half  that  of 
the  projectile  with  area  ratio  of  0.7.   The  higher  drag  co- 
efficient was  attributed  to  the  underexpanded  flow  at  the 
exit  resulting  in  high  base  drag. 

Range  tests  conducted  for  the  concept  evaluation  of 
the  20-mm  tubular  projectile  for  the  Vulcan  Weapons  System  [5] 
showed  that  the  tubular  projectile  performed  better  than  the 
standard  HEI  round,  with  lower  drag,  shorter  time  of  flight  for 
a  given  range  (30  percent  shorter  at  1000  meters,  40  percent 
shorter  at  2000  meters)  and  similar  dispersion  characteristics 
on  target. 

Recent  tests  at  the  Naval  Weapons  Center,  China  Lake, 
California  of  the  Ball  Obturated  Tubular  Projectile  (BOT) 
have  shown  a  definite  advantage  in  that  the  low  drag  of  the 
tubular  projectile  means  slower  retardation  of  velocity  and 
hence  a  higher  terminal  kinetic  energy  [6].   (The  BOT  is  the 
object  of  this  investigation  and  will  be  described  in  detail 
in  later  sections  of  this  thesis.) 


;7 


2 .   Target  Impact  Tests 

The  higher  the  energy  delivered  to  the  target,  the 
greater  the  damage.   Since  kinetic  energy  is  directly  propor- 
tional to  the  first  power  of  mass  and  to  the  square  of  the 
velocity,  the  reduced  mass  of  the  tubular  projectile  is  com- 
pensated by  a  higher  muzzle  velocity  and  therefore  a  potential 
for  higher  kinetic  energy.   However,  due  to  the  improved  drag 
characteristics,  the  reduction  in  mass  to  obtain  the  higher 
muzzle  velocities  may  not  be  necessary,  depending  upon  the 
imposed  performance  criteria. 

Target  impact  studies  have  been  carried  out  with 
various  projectile  configurations  and  target  types. 
Rhethorst,  et.al.  [7]  conducted  impact  studies  of  7.62-mm 
tubular  projectiles  on  helmets.   Tests  showed  that  even  with 
the  same  energy,  the  tubular  projectile  penetrated  further. 

Kitchen  and  Keeser  [8]  conducted  studies  for  the 
Air  Force  on  the  impact  effectiveness  of  tubular  projectiles 
on  simulated  aircraft  fuel  cells.   These  tests  were  conducted 
with  steel  and  depleted  uranium  (DU)  tubular  projectiles  and 
standard  20-mm  HEI  projectiles.   The  projectiles  were  fired 
at  double  panels  at  varying  degrees  of  obliquity.   Of  the 
forty-seven  tests  of  the  steel  tubular  projectiles  which  im- 
pacted the  target,  twenty  breached  the  rear  panel  up  to 
angles  of  70  .   The  DU  projectiles  breached  the  rear  panel 
even  at  angles  of  85°  and  fires  were  started  in  three  of  the 
five  DU  tests.   The  standard  M56  HEI  projectiles  failed  to 
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breach  the  rear  panel  in  each  of  the  five  firings  even  though 
severe  damage  to  the  front  panel  and  two  fires  occurred. 

Brunsvold  and  Kalivretenos  [9]  conducted  a  program 
to  test  the  effectiveness  of  the  tubular  projectile  against 
a  simulated  cruise  missile  warhead.   Tests  were  conducted 
with  20  and  40-mm  rounds.   Results  of  those  tests  are 
classified,  and  therefore  not  presented  here.   It  was  noted 
[9]  that  due  to  the  improved  stability  characteristics  of 
the  tubular  projectile,  only  half  the  spin  rate  need  be  im- 
parted by  the  rifling  in  the  barrel  to  obtain  the  same  sta- 
bility as  a  standard  projectile,  which  indicates  less  barrel 
wear. 

3.   Weapons  System  Capability 

For  the  tubular  projectile  to  be  practical,  it  must 
be  able  to  interface  with  existing  weapons  systems.   Refer- 
ence 5  is  a  report  of  the  19  7  8  tests  conducted  by  the  Army 
to  evaluate  the  tubular  projectile  in  the  Vulcan  weapons 
system. 

Because  the  nose  of  the  tubular  projectile  is  flat 
compared  to  the  standard  projectile  (See  Fig.  1),  it  " .  .  . 
did  not  lend  itself  to  chambering  in  the  weapon."   Personnel 
safety  required  that  the  weapon  be  remotely  operated  and 
therefore  only  the  surface-to-surface  performance  comparison 
was  made. 

The  dispersion  on  target  improved  for  the  tubular 
projectile  from  700  to  2000  meters.   The  dispersion  varied 
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Figure  1.  20-min  Ball  Obturated  Tubular  Projectile  (BOT)  , 
components,  and  conventional  20-min  projectile. 
(Photograph  courtesy  of  NWC,  China  Lake) . 


Figure  2a. 


Standard  20-iran  projectile  in  flight  with 
detached  bow  shock.   (Photograph  courtesy 
of  NWC,  China  Lake) . 
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Figure  2b. 


BOT  with  passage  blocked  resulting  in 
detached  bow  shock.   (Photograph 
courtesy  of  NWC,  China  Lake) . 


Figure  2c. 


BOT  with  no  blockage  in  passage;  only  attached 
oblique  shocks  emanating  from  the  lip  of  the 
projectile.   (Photograph  courtesy  of  NTVC,  China 
Lake) . 
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on  the  average  by  only  0.2  mil  over  the  entire  range  spectrum 
between  the  tubular  and  conventional  HEI  projectile.   However, 
at  2000  meters,  the  tubular  was  better  by  0.4  mil.   The  time 
of  flight  to  target  was  another  impressive  finding  in  Ref.  5. 
The  tubular  projectile  had  a  30  percent  shorter  flight  time 
at  1000  meters  and  a  40  percent  shorter  flight  time  at  2000 
meters  than  the  conventional  HEI  projectile. 

Results  of  firings  of  25-mm  and  30-mm  tubular  pro- 
jectiles from  the  Oerlikin  KBA  and  GAU-8/A  respectively  were 
reported  in  part  in  Ref.  3.   Figures  comparing  Line-of-Sight 
Penetration  Capability  vs.  Range  for  a  conventional  spinner 
configuration  based  on  AR-2  shape,  a  f inner  configuration 
based  on  FlO  shape,  and  a  tubular  projectile  were  shown.   The 
tubular  projectile  was  shown  to  be  inferior  in  both  instances. 
The  report  does  not  indicate  the  degree  of  compatibility  the 
tubular  configuration  has  with  the  guns  used. 

4.   Review  Of  Internal  Shock  Wave  Considerations 

From  photographs  of  tubular  projectiles  in  flight  [4] 
and  during  wind  tunnel  testing  [7],    various  shock  patterns 
have  been  observed  at  the  inlet  region  and  in  the  wake.   The 
photographs  of  Figs.  2a-c  show  the  BOT  and  the  conventional 
projectiles  in  flight.   The  sabots  used  in  Refs.  4  and  7  were 
of  the  pusher  type  similar  to  those  illustrated  in  Fig.  3. 
The  detached  bow  shock  seen  in  the  photograph  of  the  standard 
projectile  (Fig.  2a)  is  similar  to  that  of  the  tubular  pro- 
jectile shown  in  Fig.  2b  with  passage  closed  (and  in  Ref.  4 
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with  the  sabot  attached) .   With  no  blockage  of  the  internal 
passage  (Fig.  2c)  there  was  no  detached  bow  shock,  with  only- 
attached  oblique  shocks  emanating  from  the  lip  of  the  pro- 
jectile.  This  latter  case  allows  the  possibility  of  four 
other  internal  flow/shock  configurations. 

(1)  A  normal  shock  standing  at  the  entrance  to  the  pro- 
jectile; 

(2)  A  normal  shock  standing  at  some  intermediate  position 
in  the  channel; 

(3)  A  system  of  oblique  shocks  present  in  the  channel; 

(4)  The  channel  is  devoid  of  shocks  and  the  flow  is  shock- 
free  through  the  channel. 

In  each  of  the  cases  described,  specific  conditions 
must  exist: 

a.   The  Detached  Bow  Shock 

When  the  projectile  passage  is  closed,  such  as 
in  the  case  of  the  attached  sabot,  the  projectile  acts  as  a 
blunt-nosed  body  with  the  resultant  detached  normal  shock  [10] 
(see  Figs.  2a  and  b) .   From  oblique  shock  theory,  there  is  a 
maximum  turning  angle  which  the  flow  (supersonic)  can  nego- 
tiate through  an  attached  oblique  shock.   When  this  turning 
angle  is  exceeded,  the  shock  becomes  normal  and  detaches. 
There  have  been  many  studies  as  to  the  strength  and  location 
of  this  detached  bow  shock,  but  here  the  internal  flow  is  of 
more  interest  and,  therefore,  those  references  are  not  in- 
cluded in  this  discussion. 
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b.   Normal  Shock  At  The  Entrance 

From  one-dimensional  frictionless  flow  theory, 
for  a  given  initial  Mach  number,  M  ,  the  pressure  ratio 
across  a  normal  shock  is  given  by: 


P     (1  +  yM^) 

-i  =  L  (1) 

P,    (1  +  ym') 

2 


for  a  perfect  gas  with  constant  specific  heats.   P   is  the 

2 

pressure  downstream  of  the  normal  shock,  P   the  pressure  up- 
stream of  the  shock,  M  the  downstream  Mach  number  and  y  the 

2 

ratio  of  specific  heats,  C  /C^.   This  pressure  ratio  must 
exist  in  order  for  a  normal  shock  to  be  formed.   Kantrowitz 
and  Donaldson   [11]   conducted  experiments  in  which  they  con- 
cluded that  the  contraction-ratio  in  the  entrance  of  a  super- 
sonic diffuser  was  important  in  reducing  the  kinetic  energy 
losses.   The  studies  done  were  on  supersonic  diff users  with 
varying  areas.   The  formation  of  a  normal  shock  at  the  en- 
trance is  determined  by  the  throat  area  and  initial  Mach 
number  such  that  the  Mach  number  at  the  throat  of  the  duct 
is  exactly  1.0.   The  maximum  value  of  the  ratio  of  the  capture 
(inlet)  area,  A^,  to  throat  area,  A^ ,  required  to  allow  the 
bow  shock  to  be  swallowed  is  given  in  Ref.  9  as: 

At    M    L  Y+1   ^2    2^'' 
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where  M  may  be  obtained  from  tabulated  values  for  normal 

2 

shock  waves  for  a  given  free  stream  Mach  number,  M  or  may 

1 

be  calculated  from 

2 


M   = 

2 


(Y-1) 


+  M^ 


2y     5 


(3) 


L(Y-i)  '\ 

c.   Normal  Shock  At  An  Intermediate  Position 

The  formation  and  position  of  a  normal  shock  at 
an  intermediate  position  in  a  constant  area  duct  has  been 
theorized  in  the  presence  and  absence  of  viscous  effects. 
In  the  absence  of  viscous  effects,  the  normal  shock  tends 
to  grow  from  the  interaction  of  oblique  shocks  and  will  be 
discussed  under  "d"  below. 

In  the  presence  of  friction,  the  formation  of 
the  boundary  layer  along  the  duct  walls  reduces  the  effective 
area  of  the  duct  [12,  13,  14].   This  reduction  has  the  effect 
of  decelerating  a  supersonic  flow  as  in  a  diffuser.   In  a 
constant  area  duct  with  friction,  the  location  of  a  normal 
shock  in  the  duct  may  be  determined  by  Fanno  Line  Flow  Theory. 
The  length  to  diameter  ratio,  friction  factor  and  initial 
Mach  number  determine  the  range  of  back  pressures  for  which 
the  duct  will  be  "choked",  or  a  normal  shock  set  up  [10]. 

For  a  given  supersonic  entrance  Mach  nijimber,  there 
exists  a  maximum  value  of  fL/D;  f,  being  the  Reynolds  number- 
dependent  friction  factor,  L,  the  length  of  the  duct  and  D, 
the  diameter  of  the  duct.   This  maximum  value,  is  the  value 
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for  which  the  duct  will  remain  unchoked.   If  the  actual  value 
of  fL/D  is  less  than  the  maximum,  the  flow  will  remain  un- 
choked if  the  back  pressure  is  sufficiently  low  at  the  duct 
exit. 

In  the  case  of  the  tubular  projectile,  assuming 
a  constant  back  pressure,  L/D  of   (typically)  8.0  and  Re  in 
the  range  from  approximately  4x10^  to  6x10^,  from  Fig.  20.1 
of  Ref.  15,  f  is  between  0.013  and  0.0125.   Using  an  average 
value  of  f  of  0.01275, 

fL/D  =  0.102 
From  Ref.  10,  Appendix  E,  this  corresponds  to  M=:1.41. 
Iterating  to  obtain  better  values  of  f  from  a  revised  average 
Re  would  increase  the  Mach  number  at  which  the  projectile 
would  choke.   After  one  iteration  the  Mach  number  increases 
to  M  =  1.44.   As  L/D  decreases,  Re  increases  and  f  decreases 
resulting  in  a  lower  value  of  the  Mach  number  at  which  the 
flow  chokes. 

d.   System  Of  Oblique  Shocks  In  The  Duct 

The  formation  of  a  system  of  oblique  shocks  in  a 
duct  or  channel  is  the  most  difficult  to  explain  and  analyze. 
Ferri  [16,  17]  has  developed  an  analytical  solution  for  the 
shock  structure  utilizing  the  Method  of  Characteristics. 
The  analysis  becomes  difficult  if  not  impossible  at  the 
centerline  of  the  duct  because,  according  to  Ferri,  the  oblique 
shock  curves  toward  the  centerline  and  becomes  normal.   The 
extent  of  the  normal  shock  is  dependent  upon  geometry  and  flow 
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conditions.   The  larger  the  wedge  angle  for  a  given  Mach 
number,  the  stronger  and  larger  the  normal  shock  at  the 
centerline.   The  constant  area  duct  implies  a  zero  internal 
deflection  angle.   In  this  instance,  Ferri's  analysis,  which 
ignores  viscous  effects,  implies  that  no  normal  shock  is 
formed  (Fig.  13  of  Ref.  [16]).   With  no  wedge  angle,  there 
is  no  oblique  shock  formed. 

Brunsvold  and  Kalivretenos  [9]  have  modified  a 
Method  of  Characteristics  computer  code  [18]  developed  for 
supersonic  inlets  to  account  for  the  absence  of  a  center-body 
on  the  tubular  projectile.   This  modified  code  was  then 
coupled  with  a  subroutine  to  calculate  the  internal  wave  drag 
coefficient  based  on  A   and  used  to  obtain  an  optimum  internal 
geometry  for  the  tubular  projectile  which  would  give  minimum 
internal  wave  drag.   The  computer  code  assumes  inviscid, 
supersonic  flow  throughout.   Reference  9  states  that  some 
inaccuracies  are  allowed  in  the  computation  due  to  the  re- 
laxation of  the  convergence  of  the  characteristics  equations 
at  the  centerline,  but  that  the  inaccuracies  do  not  have  time 
to  propagate  to  the  walls  of  the  projectile  before  the  exit 
plane.   Therefore,  the  inaccuracies  do  not  affect  the  wall 
pressure  distribution. 

Rhethorst,  et.al.  [19]  has  postulated  another 
method  for  determining  the  shock  structure  and  the  resultant 
internal  wave  drag  for  a  tubular  projectile,  using  Prandtl- 
Meyer  wave  theory.   But  again  the  initial  oblique  shock 
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formation  is  dependent  upon  the  wedge  angle  at  the  entrance 
and  therefore  does  not  apply  to  the  constant  area  duct  case 
where  no  internal  wedge  angle  exists.   Ferri  [16,  17], 
Brunsvold  and  Kalivretenos  [9]  and  Rhethorst,  et.al.  [19] 
have  neglected  boundary  layer/shock  wave  interaction. 

Fejer,  et.al.,  [14]  and  Waltrup  and  Billig  [15] 
have  conducted  experiments  utilizing  constant  area  ducts. 
Fejer,  et.al.,  performed  experiments  on  constant  area  rec- 
tangular ducts  varying  in  lengths,  whereas  Waltrup  and  Billig 
conducted  their  experiments  on  cylindrical  ducts  of  varying 
lengths  and  diameters.   Other  than  the  obvious  geometry  dif- 
ferences, the  procedures  were  very  similar.   The  advantage 
of  the  rectangular  duct  was  the  ability  to  visualize  the  shock 
patterns  through  a  viewing  port,  whereas  the  curvature  of  the 
walls  of  the  cylindrical  duct  precluded  accurate  visualization. 

In  both  experiments,  oblique  shock  systems  were 
formed  in  the  ducts.   The  location  of  the  shock  system  in  the 
duct  and  the  actual  configuration  of  the  shock  system  was  con- 
trolled by  the  pressure  ratio.   The  pressure  ratios  were  con- 
trolled in  different  ways,  but  the  effects  were  the  same.   As 
back  pressure  increased,  a  shock  system  formed  at  the  exit  of 
the  duct.   With  further  increases  in  back  pressure  the  system 
of  shocks  moved  upstream  and  out  of  the  duct  and  a  detached 
shock  appeared.   Shapiro  [12]  indicates  that  internal 
oblique  shocks  may  be  the  remnants  of  normal  shocks  "with 
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bifurcated  ends  as  a  result  of  boundary  layer  separation." 
He  states  that  the  normal  portion  grows  shorter  as  the 
boundary  layer  gets  thicker  until  the  normal  portion  dis- 
appears completely  leaving  only  the  bifurcated  portion.   He 
illustrates  this  in  Fig.  28.27(c)  of  Ref.  12.   This  explains 
the  movement  of  the  oblique  system  toward  the  entrance  of  the 
duct  in  the  experiments  of  Refs.  14  and  15.   As  back  pressure 
increases,  the  flow  downstream  of  an  internal  shock  has  more 
of  a  tendency  to  separate.   Fejer  also  states  that  the 
presence  of  a  turbulent  boundary  layer  had  a  stabilizing 
effect  on  the  shock  system  in  the  constant  area  rectangular 
duct.   In  the  experiments  involving  the  cylindrical  duct,  the 
shock  structure  again  moved  from  exit  toward  entrance  with 
back  pressure  increases.   They  also  report  a  breakdown  in  flow 
downstream  of  the  shock  system  as  indicated  by  relatively 
constant  pressure  movements  in  this  region, 
e.   No  Shocks  Present  In  The  Duct 

It  can  be  seen  from  previous  discussion  that 
given  an  initial  Mach  number  and  inlet  stagnation  pressure, 
a  sufficiently  high  back  pressure  will  cause  a  normal  shock 
or  system  of  oblique  shocks  to  be  formed.   If  the  cross- 
sectional  area  of  the  duct  is  reduced  in  the  direction  of 
flow,  as  in  a  supersonic  diffuser,  the  incident  flow  may 
cause  an  oblique  compression  wave  to  be  formed.   In  the 
absence  of  the  requisite  back  pressure  or  area  reduction. 
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no  explanations  for  internal  shock  formation  have  been  found 
other  than  those  due  to  viscous  effects  noted  in  "d"  above. 

f .   Summary 

In  a  constant  area  duct  with  supersonic  inlet 
flow,  there  is  a  pressure  rise  caused  by  a  reduction  in  the 
effective  flow  area  by  boundary  layer  formation  due  to 
friction.   Depending  upon  the  length-to-diameter  ratio  (L/D) 
of  the  duct  and  the  flow  parameters,  the  increase  in  pressure 
ratio  (or  reduction  in  Mach  number)  required  to  cause  forma- 
tion of  a  system  of  oblique  shocks  is  less  than  that  required 
to  set  up  a  normal  shock  at  the  duct  entrance.   This  difference 
is  explained  by  the  boundary  layer  formation  and/or  interaction 
with  the  shock  waves. 

In  the  case  of  the  tubular  projectile,  the  back 
pressure  at  the  exit  in  flight  is  affected  by  the  complex  wake 
flow.   There  is  no  pressure  difference  in  the  surrounding 
medium  at  large  distances  from  the  projectile.   The  formation 
of  shock  waves  inside  the  projectile  with  a  constant  cross 
section  is  due  to  viscous  effects  and  depends  upon  Mach  number 
and  the  length-to-diameter  ratio. 

5.   Laboratory /Wind  Tunnel  Testing 

From  all  indications,  the  Weapons  System  Concept  Team 
(WSCT)  of  ARRADCOM  has  done  extensive  wind  tunnel 
testing  on  various  tubular  projectile  configurations.   The 
work  of  A.  Plateau  of  the  Weapons  System  Concept  Office  is 
referenced  often  with  regard  to  results  obtained,  however  no 
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publications  from  that  office  have  been  found.   (The  work 
referenced  in  other  publications  is  usually  of  the  "private 
communication"  type.)   Appendix  D  of  Ref.  8  reports  results 
of  a  test  conducted  by  the  WSCT  at  the  ARRADCOM  facility. 
The  facility  utilizes  a  variable  Mach  number  tunnel  which 
may  be  used  to  simulate  the  deceleration  of  the  projectile 
in  flight. 

The  results  of  Ref.  8  agree   closely  to  those  of 
Ref.  1 ,    in  that  the  drag  coefficient  rises  almost  instan- 
taneously at  the  Mach  niomber  at  which  the  flow  becomes 
choked. 

The  Mach  number  at  which  this  choking  occurs  is  de- 
pendent upon  geometry.   For  the  case  of  inviscid  flow, 
Brunsvold  and  Kalivretenos   [9]   give  the  following  relation 
for  the  area  ratio  corresponding  to  the  free  stream  Mach 
number  at  which  the  shock  is  "regurgitated"  and  the  flow  be- 
comes choked: 


y+1 
^2(Y-1) 

—  =   -  I  A  f  1  +  ^  M''  ^ 
At       Ml   Y+1   I      2    1 


(4) 


C.   PURPOSE  OF  STUDY 

The  results  of  tests  and  experiments  previously  mentioned 
show  for  the  most  part  definite  advantages  to  the  tubular 
projectile  over  the  conventional  projectile  primarily  in  the 
area  of  lower  drag  (approximately  1/2  to  1/3  of  conventional  [s]) 
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which  allows  a  flatter  trajectory  and  shorter  time  of  flight 
for  a  given  range  [5 J .   With  the  emphasis  on  Anti-Missile 
Defense  both  from  a  Surface-to-Air  and  Air-to-Air  standpoint, 
these  features  are  definitely  worth  further  study. 

The  problem  seen  from  the  air  platform  is  that  of  the 
discarding  sabot  or  pusher.   The  solid  [4,  5,  7]  or  split  [8] 
disk  is  effective  and  reliable,  however  the  possibility  of 
ingestion  in  the  aircraft  engines  prohibits  their  use. 
Rhethorst  et.al.  in  Ref.  7  illustrates  a  number  of  sabot/ 
obturator  designs  for  use  with  the  tubular  projectile,  but 
these  are  also  of  the  discarding  type.   They  also  investigated 
the  possibility  of  using  a  consumable  sabot  [19]  which  would 
be  burned  up  as  the  projectile  was  transiting  the  barrel  and 
would  be  completely  consumed  as  it  exited  the  barrel.   The 
most  recent  development  has  been  the  BOT  which  was  designed 
at  NWC  China  Lake,  California  [6].   The  ball  obturator  has 
been  bored  with  a  hole  the  same  diameter  as  the  hole  through 
the  projectile.   When  loaded,  the  ball  is  supported  inside 
the  projectile  such  that  the  axis  through  the  hole  in  the  ball 
makes  an  angle  of  approximately  90°  with  the  axis  of  the  hole 
through  the  projectile  (see  Fig.  4) . 

Gas  pressure  from  the  burning  propellant  holds  the  ball 
obturator  fixed  with  respect  to  the  projectile  as  the  rifling 
in  the  barrel  imparts  spin  to  the  projectile  (as  well  as  the 
ball).   Upon  exiting  the  barrel,  the  propellant  gas  pressure 
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is  released.   At  that  time  a  complex  force  distribution  acts 
on  the  ball  to  cause  it  to  change  position  inside  the  pro- 
jectile in  such  a  way  as  to  align  the  holes.   Inertial  forces 
then  dominate  to  maintain  this  alignment  so  that  the  pro- 
jectile remains  fully  tubular.   There  are  no  separating  parts 
and  the  opening  process  is  automatic. 

The  purpose  of  this  research  is  to  analytically  predict 
the  motion  of  the  ball  obturator  as  a  function  of  time.   In 
support  of  this  purpose,  an  experimental  apparatus  was  de- 
signed and  built  to  simulate  the  motion  of  the  obturator  in 
order  to  gain  physical  evidence  of  the  accuracy  of  the  model. 
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II.  ANALYSIS 

A.   DEVELOPMENT  AND  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

The  ball  obturator  may  be  thought  of  as  a  rigid  body  of 
rotation  with  a  system  of  coordinate  axes  fixed  to  the  ball 
having  its  origin  at  the  mass  center  of  the  ball  (Fig.  5). 
These  axes  will  be  designated  x,  y,    and  z,    where  z  is  the 
axis  through  the  hole  in  the  ball. 

The  motion  of  the  ball  jnust  be  described  relative  to  the 
projectile,  and  for  this  purpose  an  inertial  frame  of  refer- 
ence is  defined  with  axes  designated  X,  Y  and  Z  with  its 
origin  also  at  the  mass  center  of  the  ball.   For  this  analysis 
it  is  assumed  that  the  projectile  has  only  a  rotation  about 
the  Z-axis  and  no  translation  or  other  rotation.   Because  only 
the  motion  of  the  ball  with  respect  to  the  projectile  is  de- 
sired at  this  point,  the  relative  velocities  are  the  quantities 
desired.   The  ball  has  only  rotational  components  and  thus  the 
values  of  angular  velocity  with  respect  to  the  X  and  Y  axes 
are  also  relative  to  the  projectile,  whereas  the  value  with 
respect  to  the  Z-axis  must  account  for  the  projectile  spin. 

By  doing  this,  the  motion  may  be  described  by  Euler's 
Modified  Equations  of  Motion  for  a  rigid  body  of  revolution 
about  a  fixed  point.   These  are  given  in  Thomson  [20]  as: 
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Aco   +  (C-A)w  0)   =  G  (5a) 

12   3       1 

Ad)   +  (A-C)a)  to   =  G  (5b) 

2  13       2 

Co)   =  G  (5c) 

3       3 


where 

A        =  Mass  moments  of  inertia  about  principal 
axes  perpendicular  to  the  z-axis. 

C        =  Mass  moment  of  inertia  about  the  prin- 
cipal axis  through  the  hole  (z-axis) . 

Ci    ,(h    ,11^     =  Angular  accelerations  about  the  body 
^   ^   ^    fixed  axes  x,    y,    z  respectively. 

03  ,0)  ,0)   =  Angular  velocities  about  the  body  fixed 

12    3 

axes. 
G  ,G  ,G  =  Moments  acting  about  the  body  fixed  axes. 

12    3 

Since  a  body  of  revolution  has  only  two  distinct  prin- 
cipal axes,  the  mass  moments  of  inertia  about  the  x  and  y 
axes  are  equal.   The  moments  of  inertia  of  the  ball  were 
determined  by  first  finding  the  inertia  tensor.   From 
Refs.  21  and  22  the  inertia  tensor  is  found  by  starting  with 
the  formulation  for  a  hollow  cylinder  of  inside  radius  r  and 
length  2h  (Fig.  6)  in  which  the  inertia  tensor  in  terms  of 
dyadics  is 

I  =  (U  +  ee)mrV2  +  (U-ee)mh^/3  (6) 

In  Eq.  (6)  m  is  the  mass  of  the  cylinder  (ball) ,  U  is  the 
unit  dyadic  and  ee_  is  the  referenced  dyad.   Using  Fig.  6  as  a 
guide,  the  incremental  value  of  the  inertia  tensor  with  respect 
to  the  center  of  mass,  dl^  is  given  by: 

38 


(  h  =  /R^-p2  ) 


Figure  6.   Dimensional  relationships  for  calculating  the 
inertia  tensor  and  moments  of  inertia. 
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dig  =  (U  +  ee)pV2  dm  +  (U  -  ee)hV3  dm  (7) 

where  p  is  the  radial  distance  to  the  incremental  mass,  dm. 
The  incremental  mass  is  given  by 

h 
dm  =  47Tp(R2-p2)  Ydp  (8) 

From  Eq.     (7) ,    since   the   origin   and   center  of  mass    coincide, 

I^  =    /   dl^  (9a) 

o        J         G 

therefore , 

•R  rR 


I      =   /     (U+ee)  2p^TThYdp    +   /      (U-ee_)    2—  47TpYdp 
-'r  -^r 


which,    after   simplification    leaves 


I      =    27TY(U+ee)         p^(R^-p^)    dp   + 


•R  3/2 

3 -'' 


4-7TY(U-ee)/     (R^-pM  pdp 


(9b) 


(9c) 


Solving  Eq.     (9c)    gives 

3/2 

I      =   7rY{2(U+ee)  [|:r^+|-R^(R^-r2)         ]    + 
o  —      b         15 
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Eq.  (9d)  is  the  expression  for  the  second  order  inertia 
tensor  of  the  ball  obturator  with  outside  radius  R,  hole 
radius  r  and  made  of  a  material  of  density  y.      Diadics  have 
been  used  to  describe  this  tensor.   In  order  to  obtain  the 
moments  of  inertia  about  any  axis,  the  diadic  describing  the 
tensor  need  only  be  premultiplied  and  postmultiplied  by  scalar 
multiplication  by  the  unit  vector  along  the  axis  about  which 
the  moment  of  inertia  is  described.   It  may  be  seen  in  Fig.  6 
that  the  e-direction  vector  is  along  the  z-axis.   The 
e-direction  may  be  chosen  arbitrarily,  however  this  choice 
greatly  simplifies  the  formulation  of  the  integral.   Thus 
the  moment  of  inertia  about  the  z-axis,  C,  is  given  by 


1    2  3/2 

C  =  e  •  [  (ii+jj-t-ee)+ee)  ]  *  e  {  [2117  ( jr^+j^-R^  (R^-r^)    ] 


4  •     ;L       5/2 

+  3(i£+2l+ee-ee)  [iryjCR^-r^ )    ]}  (10a) 


Simplifying  ^nd  using  the  subscript  s  to  denote  particular 

values  for  the  ball  obturator: 

3/2 

C  =  4-Kyi-^Zs^    +   Xs-Rg')  (Rs'-rs')  (10b) 


For  the  axes  perpenducular  to  e  through  the  mass  center,  the 
Moment  of  Inertia  is  given  by 


A  =  B  =  j-I^-j  =  i-I^-i  (11a) 
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Performing   steps    as    in  Eqs.     (10)    gives 


A=    j-  (ii+j  j+ee+ee)  »  j     [2irY  (ir^+i^R^  (R^-r^  )    '     ] 


.4 


3/2 

5/2 


+   j-  (ii+Ji+ee-ee)  •  j     [TTY-^CR^-r^)         ]  (lib) 


After  premultiplying  and  postmultiplying. 


1      .         2         _       ,         .    ^/2.         4  .         .5/2 

' (lie) 


A=    TTY{2[5rs^    +   JS^^s^  (Rs^-^s^)  ]    +   l5(Rs^"^s^^  ^ 


In  order  to  simplify  Eqs.  (5),  define  a  non-dimensional 
X  such  that 

X  =  ^  (12) 


Thus   Eqs.     (5)    becomes 

d)      +  Xo)   CO  =   G  /A  (5a) 

1  2      3  1 

0)      +   Xco   0)  =    G   /A  (5b) 

2  13  2 

(L  =    G   /C  (5c) 

3  3 


The  position  of  the  ball  at  any  time  may  be  described  by 
three  angles,  0,  (^   and  \\>,      These  are  known  as  the  Euler  angles 
and  their  relationships  to  the  coordinate  axes  already  de- 
scribed may  be  seen  in  Fig.  7.   Reference  20  defines  oj  ,  o) 
and  to   in  terms  of  the  Euler  angles  as 

3 
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0)  =  i);sinesin4)   +  5cos(J)                                                                        (13a) 

oj  =  ii;sinecos(i)   -  esin(f)                                                                        (13b) 

(0  =  <})   +   ipcosG  (13c) 

3 


By  differentiating  Eqs.     (13)    with   respect   to   time,    iji    ,    di    , 

1  2 

and  (j      may  be   obtained  as 

3 


0)      =   ipsinesmcj)   +   i|;  ( GcosSsincf)   +    (|)sinecos(}))    +    BcoscJ)         (14a) 
-e^sincj) 

(I)      =   'I^'sin9coS(|)    +  tJ  (9cos9cos4>    -   isinSsintj) )    -   9sin(})         (14b) 

-    9^cos({) 
^      =    $   +   i|;cos9    -   i|^9sin9  (14c) 

3 

By  substituting  Eqs.  (13)  and  Eqs.  (14)  into  Eqs.  (5), 
the  following  expressions  for  G  /A,    G  /A  and  G  /C  are  found 

12  3 

to  be 


i|!isin9sin(J)   +   i|;  ( 9cos9sin(})   +    ^isinGcoscj))    +    Scoscf)    -  (15a) 

•    •  •  •  •  • 

9(|)sin(l)   +    X  (4)sin9cos4)    -    9sin(j))  (\jJCOs9    +    <p)    =   G   /A 

i|;sin9cos(J)   +   i];  ( 9cos9cos(})   -    4>sin9sin4))    -    9sin({)   -  (15b) 

94)cos(J)   -    X  (ij;sin9sin(|)  +    9cos(J))  (ij;cos9   +    i)    =   G  /A 

ipcosG   -   ijJG'sinG   +    $  =   G  /C  (15c) 

3 

Transforming  this  system  of  differential  equations  into  the 
matrix  form   [a]  {x}  =  [B]  leaves 
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sinesin(J)      0        coscj) 

sinecos(})      0      -sin(|) 

cos 9  1  0 


$      = 


G  /A   -    ij;  Ce  cose  sin  (|)    +    ^sinecos(|)) 

•         •  • 

G   /A   -    1^  (e  cose  cos  (|)    -    ({)sinesin4)) 

2 

G  /C   +   ij;esine 

3 


+    e^Jsin^   -    X  (ii;sinecos(|)   -    esincji)  (^cose   +    (^) 
+   ej)cos4)  +   X  (lysine sin (|)  +   ecos(j))  (^cose   +   J) 


(16) 


By  premultiplying  both  sides  of  Eq.  (16)  by  the  inverse 
of  the  coefficient  matrix,  we  are  left  with 


G  /A  -  i|;  ( 5cosesin(j)  +  ^sinecoscj)) 
1 


sin  4) 

cosip 
sine         ^ 

sine 

$ 

>  = 

sin<j)cose 
sine 

-cosecos(t) 
sine         ■'• 

\    J 

cos^ 

-sincj)         0 

G  /A  -  4^  (e  cose  cose})  -  (|)sinesin(})  ) 

2 


G  /C   +   ijjesincj) 

3 


+    e^sinci)   -    X  (ijjsinecosct)   -    esin(|))  (ij^cose    +    ^) 
+    9i)cos(j)   +    X  (ij;sinesin({)   +    ecosc}))  (ijjcose   +    }>) 


(17) 


In  order  to  solve  this  system  of  non-linear  second  order, 
differential  equations  for  e,  (j)  and  ^,    the  following  must  be 
known : 

Initial  Values  (denoted  by  subscript  o)  of  e,  ^,    ^,    Q,    <p 
and  i  and  the  quantities.  A,  C,  G  ,  G   and  G  . 
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Due  to  symmetry,  0°  £  0  £  90°,  and  since  we  are  interested 
in  the  performance  of  the  ball  away  from  0°,  9   >  0.   Since 
the  ball  is  attached  to  the  projectile  initially  (no  slippage) , 

•         •  • 

6  =  (i)Q  =  0 .   ipQ  ■'■^  equal  to  the  spin  rate  of  the  projectile 
as  it  exits  the  barrel.   A  and  C  are  physical  constants  de- 
pendent upon  geometry  and  material  properties  [Eqs.  (10)  and 
(11)  ]  and  therefore  X  is  a  constant.   <J>q,  and  ipQ  are  arbitrary 
due  to  symmetry  and  the  values  of  G  ,  G   and  G  must  be  de- 

1      2  3 

termined  from  knowledge  of  forces  acting  on  the  ball. 

The   solution  to  the  problem  of  the  motion  of  an  axi- 
symmetric  rigid  body  about  a  fixed  point  with  no  moments 
applied  is  well  known  [20,21,22,2  3].   Under  these  conditions 
there  is  no  nutation  and  the  precession  rate  ^  is  constant. 
Thus 

9  =  0 

and 

i  =  ^ (18) 

^    (A-C)cos9 

The  computer  solution  method  was  tested  by  substituting  the 

projectile  spin  velocity  for  ^   in  Eq.  (18)  and  calculating 

o 
(^   for  the  specified  value  of  0  <  9q  <  90  .   By  inserting 

these  values  of  9^,  ^   and  '^   into  the  computer  model,  and 

solving  for  9  as  a  function  of  time,  9  was  indeed  found  to 

be  zero.   In  addition,  all  other  criteria  for  the  test  case, 

i.e.,  constant  angular  momentum  and  constant  ^   and  cj)  were  met. 
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The  problem  remaining  is  one  of  modeling  the  forces  acting 
on  the  ball  in  order  to  determine  G  ,  G   and  G  . 

12  3 

B.   MODELING  OF  THE  APPLIED  TORQUES 

1.   Moments  Due  To  Viscous  Effects 


Without  the  influence  of  applied  forces  or  moments, 
the  angular  momentum  remains  constant  or,  from  first  prin- 
cipals : 

r 

M  =  h  =  0  (19) 

In  order  for  the  ball  to  nutate  to  the  position  in 
which  the  holes  are  aligned,  a  moment  must  be  applied.   This 
resultant  torque  is  due  to  both  viscous  and  sliding  frictional 
forces. 

The  first  attempt  was  to  model  the  torque  on  the  ball 
as  that  due  to  laminar  Couette  flow  of  a  viscous  fluid  between 
two  flat  plates  where  the  shear  stress  is  given  by 

T   =  y.^  (20) 

The  solution  of  the  Navier-Stokes  equations  for  this 
flow  is  well  known  and  given  by 


f  =  y,  ^  (21) 

dX  f  ^y2 


where  dP/dx  is  the  pressure  gradient  in  the  direction  of 
flow  and,  for  our  calculations,  is  assumed  to  be  zero.   The 
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governing  partial  differential  equation  is  then 

=   0  (22) 


d^u  _ 


dy^ 

By  considering  the  relative  velocity,  one  plate  is  held 
fixed  while  the  other  has  velocity  of  magnitude  oo  r.   coj- 
is  the  relative  angular  velocity  about  the  instantaneous 
axis  of  rotation  and  r  is  the  perpendicular  distance  from 
the  instantaneous  axis  to  the  point  at  which  the  velocity 
is  desired. 

The  boundary  conditions  are  therefore  given  as 

u  =  -oj^r   at  y  =  0 

(23) 
u  =  0      at  y  =  h 

where  y=0  is  the  surface  of  the  ball  and  h  is  the  gap  width 
between  ball  and  projectile.   Integrating  Eq.  (22)  twice  and 
applying  the  boundary  conditions  Eq.  (2  3)  gives 

u(y)  =  -a)^r  (1-y/h)  (24) 

Thus  from  Eqs .  (20)  and  (24),  the  viscous  shear 
stress  at  the  surface  of  the  ball,  or  at  y=0  becomes 

0)  r 
To  =  ^f-^  (25) 
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The  force  due  to  viscous  shear  on  the  ball  is  given 
by  the  integral  over  the  surface  area  of  the  ball  of  the 
shear  stress  Eq.  (25) ,  or 

The  differential  area  dA  is  given  by 

dA  =  27tR_  (R.sinoda)  (27) 

which  is  shown  in  Fig.  8* 

Substituting  Eqs.  (25)  and  (27)  into  Eq.  (26)  gives 

Ftt  =     ^ /   rsinada  (28) 

^       h       j 

o 


The  value  of  r  changes  from  0  at  the  axis  of  rotation 
to  Rg  at  a=90°.   Thus 

r  =   RgSina  (29) 

Substituting  this  expression  into  Eq.  (28) 

IT 

27111  R  ^co       r 

^V  ^  f   s      r         /    sin^ada  (30) 


The  torque  imparted  to  the  ball  is  the  force  given  by 
Eq.  (30)  multiplied  by  the  moment  arm  at  which  it  acts. 
Again,  the  moment  arm  is  a  function  of  Rg  and  a   and  is  given 
by  RgSina.   Substituting  this  into  Eq.  (30)  gives  an  ex- 
pression for  the  torque  acting  on  the  ball  obturator   due  to 

c49 


Figure  8.   Differential  area,  dA,  orientation  with  respect 
to  the  relative  angular  velocity  vector,  o)  . 
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the  velocity  of  the   ball   relative   to   the   projectile    in   a 
viscous   medium,    or 

27ry   R     0)         r 
Q  =  1-^ — ^  /sin^ada  (31) 

o 

Carrying  out  this  integration  gives 

Q  =  I     ^      r   t  (32) 

3      h 

For  a  given  geometry  and  constant  fluid  properties, 
CO   is  the  only  variable  in  Eq.  (32),  therefore  that  portion 
of  the  expression  which  is  multiplied  by  oj  will  be  called 
the  shear  constant,  S,  such  that 

Q  =  S^^  (33) 


where 


S  =  |  _l_p^  (34) 


From  the  rigid  body  motion  analysis,  the  relative 
angular  velocity  is  easily  determined.   It  is 


h 

03^=    [(00  -03^  )'  +  (CO  -%  )'  +  (03  -03^  )M  (35) 

r       iPi       2P2        3P3 

where  03^  ,  a3„   and  oo^   are  the  components  of  the  projectile 
Pi    P2      P3 

angular  velocity  about  the  axes  fixed  in  the  ball.   These 
relationships  are  given  in  Thomson  [20] 
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(u   =  CO  sind)sin9  (36a) 

Pi     P 

03   =0)  cosisine  (36b) 

P2      P 

oj   =0)  cos 9  (36c) 

p3     p 

Equations  (36)  are  simplified  by  the  fact  that  a)p 
has  only  a  component  in  the  Z-direction  of  the  inertial 
frame  of  reference. 

The  moments  G  ,  G   and  G   are  obtained  by  applying 

12  3 

Eq.  (33)  and  resolving  the  moment  vector  into  components 
about  the  x,  y  and  z-axes.   The  motion  of  the  ball  may  then 
be  described  fully  by  Eq.  (17). 

Figure  9   shows  solutions  for  9  as  a  function  of  time, 
using  the  Couette  flow  model  with  the  following  inputs : 

Rg  =7.62  mm  (0.3  in.) 

rg  =  4.7625  mm  (0.1875  in.) 

h  =  0.0254  mm  (0.001  in.) 

p  =  7800   Kg/m^  (0.282  Ibm/in^  (steel) 

^air  ^  ^-^^  ^°~^  Pa-sec  (4xl0"^lbf-sec/ft2 ) 

0)  =  12042.77  rad/sec  (1. 15x10  ^RPM) 

It  appears  that  viscous  effects  alone  are  not  suffi- 
cient to  cause  the  ball  to  nutate  toward  9=0  in  a  rapid 
fashion.   To  check  the  model  again,  progressively  larger 
orders  of  magnitude  of  S  were  input  into  the  model.   The 
effect,  as  expected,  was  to  cause  0  to  approach  zero  more 
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rapidly  Csee  Fig.  10).   From  firings  conducted  at  NWC,  China 
Lake,  on  an  indoor  range,  it  is  known  that  9  reduces  to  a 
value  which  allows  light  from  behind  the  projectile  to  be 
seen  through  the  projectile  by  a  camera  aiming  at  the  pro- 
jectile image  in  a  mirror  down  range.   The  length  of  the 
range  is  such  that,  for  this  to  occur,  9  must  reach  the  value 
3g  prior  to  impact  on  target,  or  in  about  0.1  seconds.   (see 
Fig. 11  for  3g) •   The  viscous  model  described  above  does  not 
produce  the  moments  required  to  cause  sufficient  reduction  of 
9  within  this  time  of  flight. 

The  viscous  torques  on  the  obturator  were  calculated 
above  assuming  a  solid  sphere.   The  actual  obturator  has  an 
outside  surface  area  less  than  that  of  a  solid  sphere  of  the 
same  radius.   The  surface  under  the  influence  of  viscous 
shearing  stress  is  a  very  complicated  function  dependent  upon 
the  angle  0  as  shown  by  the  cross-hatched  portion  of  Fig.  H. 
Calculations  of  the  resulting  moment  are  further  complicated 
by  the  relative  effect  of  the  missing  areas  (due  to  the 
presence  of  holes)  and  their  positions  relative  to  the 
instantaneous  axis  of  rotation.   The  time-dependent  angles 
between  the  instantaneous  axis  of  rotation  and  the  z  and 
Z-axes  must  be  known  in  order  to  determine  the  limits  of 
integration  for  correct  moment  arm  consideration. 

Because  the  moment  due  to  viscous  effects  is  dependent 
upon  the  surface  area  involved  and  because. the  surface  area 
in  the  first  calculation  (assuming  a  solid  sphere)  is  larger, 
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it  is  to  be  expected  that  by  taking  the  holes  into  account 
the  true  moment  due  to  viscous  effects  will  be  reduced.   The 
effects  of  these  area  variations  will  be  "blurred"  by  the 
high-speed  relative  motion  and  the  additional  complications 
necessary  to  account  for  them  are  not  warranted  in  light  of 
the  goals  of  this  analysis.   Therefore  for  further  develop- 
ment of  viscous  effects,  the  surface  area  of  the  obturator 
will  be  considered  that  of  a  solid  sphere. 

Nakabayashi  [24]  has  conducted  experiments  to  de- 
termine the  viscous  torque  on  rotating  concentric  spheres 
with  various  fluids  and  gap  widths.   He  gives  the  general 
expression  for  the  frictional  moment  on  a  sphere  due  to  flow 
in  the  gap  as 

where  p  is  the  density  of  the  fluid  in  the  gap.   For  laminar 
flow 

C  =  TIT  — - —  Re  ^  (38) 

m   -^    ^ 


and   for   the   turbulent   flow 

Re  "0.25 
P      _    0.053    l+C7/4)£  ,-.Q, 


where  £    =  h/Rg . 

The  Reynolds  number  is  defined  as 

R  ^0) 
Re  =  -i 1  (40) 
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in  which  v  is  the  kinematic  viscosity  of  the  fluid  in  the 
gap.   For  the  narrow  gap  width  assumed  for  the  BOT,  the 
transition  Reynolds  number  may  be  approximated  by 

Re^  =  70e""^-^  (41) 

(from  Fig.  2a  of  Ref.  24.)   in  the  case  of  the  BOT,  the 
relative  angular  velocities  are  insufficient  to  assume  the 
existence  of  turbulent  flow  in  the  gap.   Thus  the  frictional 
moment  due  to  viscous  flow  in  the  gap  is  given  by 


•V  =  -  I'  -^  ^fRs'"r  ^  (^2) 


r 


For  small  values  of  e  (as  is  the  case  here) ,  Eqs .  (42)  and 
(33)  are  identical. 

2.   Moments  Due  To  Sliding  Friction 

It  has  been  noted  that  the  extreme  propellant  gas 
pressures  inside  the  barrel  are  assumed  sufficient  to  hold 
the  obturator  fixed  to  the  projectile  during  launch.   How- 
ever upon  exiting  the  barrel,  the  pressure  distribution  on 
the  projectile  is  due  to  the  aerodynamic  characteristics  of 
the  projectile.   Assuming  the  obturator  is  still  blocking 
the  projectile   duct,  a  detached  bow  shock  is  set  up  (see 
Fig.  2b).   The  pressure,  P  ,  downstream  of  this  shock  is  given 

2 

by  one-dimensional  shock  theory  from  Eq.  (1) .  It  is  assumed 
that  P  is  sufficiently  large  to  cause  the  obturator  to  con- 
tact the  rear  of  the  projectile  socket  housing  the  obturator. 
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The  forces  involved  depend  upon  the  projected  area 
of  the  hole  through  the  projectile  and  the  contact  point (s) 
of  the  obturator  and  projectile.   The  force  acting  on  the 
obturator  is  assumed  to  be  given  by 

F  =  TTr^  AP     [see  Eqs .  (1)  and  (3)]  (43) 

where  r  is  the  radius  of  the  hole  through  the  projectile 


and 


AP  =  P   -  P  (44) 

1       2 


From  sliding  friction  theory,  the  contact  area  is 
not  a  consideration;  only  the  magnitude  of  the  force  normal 
to  the  contact  surface  and  the  coefficient  of  friction  of 
the  materials  in  contact  are  important.   In  the  case  of  the 
obturator,  the  pressure  force,  F   acts  along  the  Z-axis. 
The  reaction  force  may  be  broken  into  a  component  tangent 
to  the  obturator  and  a  normal  component  acting  through  the 
center  of  the  obturator.   The  normal  component,  F  ,  multi- 
plied by  the  coefficient  of  friction  y^ ,  gives  the  sliding 
friction  force  acting  to  oppose  the  motion  of  the  obturator: 

Fg  =  Fj^Pg  (4  5a) 

where 

F^   =  FpCOsBg  (45b) 


so  that 


F^  =  F„cos3^y^  (-46) 

s    p    s  s 
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This  assumes  the  holes  through  the  obturator  and  projectile 
have  the  same  diameter  and  that  the  point  or  points  of  contact 
are  on  the  locus  of  points  described  by  the  intersection  of 
the  hole  through  the  projectile  and  the  spherical  socket 
housing  the  obturator  in  the  projectile.   This  locus  of  points 
is  a  circle  of  radius  rp  and  each  point  is  a  distance  Rg  from 
the  center  of  the  obturator.   Gravity  forces  are  neglected 
since  they  are  relatively  small  and  act  both  on  the  obturator 
and  the  projectile. 

For  a  "perfect"  fit  between  ball  and  projectile, 
the  sliding  friction  force  is  distributed  over  an  infinite 
number  of  points  and,  in  the  limit,  the  force  at  each  point 
vanishes.   This  situation  has  no  significance  in  reality  and, 
therefore,  engineering  judgement  dictates  making  some  assump- 
tions with  regard  to  the  point (s)  of  application  of  the 
friction  force  F^ . 

For  the  general  case,  a  point  of  contact  P  may  be 
described  relative  to  the  inertial  reference  frame  as  being 
at  an  angle  3„  from  the  inertial  Z-axis,  a  distance  R_  from 
the  center  of  the  obturator  and  at  an  angle  y  from  the 
X-Z  plane  measured  toward  the  positive  Y-axis  as  shown  in 
Fig.  12 .   For  this  approximation,  assume  two  contact  points 

in  the  X-Z  plane,  P^  and  P^.   P^  and  P^  are  180   apart  on 
the  circle  of  points  so  that  y=0  and  3=3^. 


59 


t 


N 


^\' 


J-l 

o 

•H 

73 

e 

C 

td 

to 

c; 

>i 

^ 

t3 

u 

0 

u 

s 

<D 

u 

(d 

(d 

0 

4-> 

+J 

c 

0) 

(U 

6 

3 

0 

-^3 

6 

Ul 

cn 

0) 

C 

:3 

•H 

cr 

C 

M 

■H 

o 

g 

4J 

^ 

(D 

c 

4-) 

0 

OJ 

•H 

'd 

■P 

^ 

0 

S-l 

•H 

0 

^^ 

M-l 

1+^ 

cn 

Cn 

a 

C 

•H 

•H 

^ 

^ 

CO 

•H 

c 

iH 

0 

Cfl 

•H 

+J 

cr> 

fd 

C 

iH 

•H 

0 

+J 

u 

fO 

iH 

u 

D 

•H 

o 

:  ^ 

rH       . 

+j 

td    W 

(U 

O    (U 

e 

u 

^1  o 

u  u 

^1  (1) 

0    0 

O      M-l    M-l 


fN 


U 
0^ 


60 


From  the  vector  analysis,  the  projection  of  one  vector 
onto  another  is  given  by  the  scalar  or  dot  product.   The 
vectors  P   and  P   are  given  by : 

1         2 

P   =  R^  (singi  +  cosBk) 

1  s      —       — 

(47) 
P  =  R   (-sin3i  +  cosgk) 

2  S  —  — 


Their  projections  onto  the  relative  angular  velocity 
vector  (jOj./  with  components  oa^^,  ojy  and  002/  are 


P  'co^  =  R<5  (co„sin6  +  oj^cosB) 


P  '0)^  =  R^  (-a3,sin3  +  co^cosS) 
2   r     s     X         z 


The  scalar  dot  product  is  also  defined  as 


(48) 


P  -(jo^  =  P  oj^cosn  (49) 

1   ^     1  ^ 


where  n  is  the  angle  between  the  vectors  P   and  cOj-.   Therefore 


from  Eqs.  (4  8)  and  (49) 


R  (o)  sin3  +  w  cosB)  =  P  oa  cosn 
s  X        z         1  r 


(50) 


R^(-tjo,sin3  +  oo^cosB)  =  P  co^cos^ 

S     X  Z  2   ^ 

The  magnitude  of  the  vectors  P   and  P   is  R  ,  so  solving 

1       2     s 


for  n  and  Q,   gives 
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n  =  cos"-^   —  sin3  +  —  cosB 

00  w 

r  r 


OJ  CO 

^  =  cos"^    -  —  sin3  +  —  cosB 

00  oo 

r  r 


(51) 


The  moment  arms  r   and  r   from  Fig.  12  are  therefore  given  by 


oo         00 
r  =  R^sin  [cos"M^^inB  +  -^os6)  ] 

1      *  00^  OOj. 

CO        00        h  (52) 

r  =  R_  [1  -  (-^ing  +  -^osB)^] 

1       S  00^  00^ 


00..  00. 

.  .5 in  rnos"^  (- 
2      S 


r  =  R^sin  [cos"^  (-  rr-sinS  +  j^cos3)  ] 


(53) 


00„  00^  ^ 

r  =  R_  [1  -  (^osS  -  -fsinS)'] 

2       *  00.^  C0„ 


The  moment  acting  on  the  obturator  due  to  the  sliding 

friction  due  to  aerodynamic  forces  are  given  by 

-► 

M   =  (r  +r  )  _£  cos3y-  ^  (54) 

sa     12   2       s  00 

Substituting  from  Eqs.  (52)  and  (53)  into  Eq.  (54)  leaves 


M„,  =  ^R^FpCOs6y^[l  -  (-^)2sin2S  -  (_z.)  ^cos^B  ]  ^  — 

Sa      2  ^  -        2         W^  00^  bi^ 


(55) 


Jo2 


It  should  be  noted  that  if  there  were  n  contact  points  that 
the  normal  force  at  each  point  would  be  given  by  F  /n. 
Equation  C54)  indicates  the  case  of  two  points  of  contact 
discussed  here. 

As  yet,  there  has  been  no  attempt  in  the  laboratory 
to  produce  the  aerodynamic  forces  present  on  the  obturator 
for  an  actual  projectile  in  flight.   The  sliding  friction 
forces  present  on  the  obturator  in  the  experiments  (described 
in  section  III)  are  due  to  gravity  alone.   The  obturator  is 
supported  by  the  projectile  simulator  (see  Fig.  24).   Because 
the  obturator  is  smaller  than  the  housing  in  the  projectile, 
a  point  contact  at  the  bottom  of  the  housing  is  assumed. 
This  point  has  inertial  coordinates  of  (-'Rq,    0,    0).   The 
vector  from  the  origin  to  the  point  is  given  by 

Pq  =  -^si  (56) 

Using  a  development  similar  to  that  for  obtaining  the  moment 
arm  for  the  friction  forces  due  to  aerodynamic  considerations, 
the  forces  on  the  model  due  to  gravity  may  be  determined. 
Again,  the  dot  product  is  given  by 

Pq-oJj.  =  -Rs^x  (57) 

which  is  the  same  as  P^co^cos^  • 

kj   r 

Therefore  the  angle  (E,)    between  the  angular  velocity 

vector  (jOj.  and  Pq  is  obtained  from 


.6  3 


C  =  cos-^  (-  -ii)  (58) 

r 


The  moment  arm  r   is  found  as  in  Eq.  (55)  from 


r^  =  R  sin  [cos'^  (-  ^)  ]  (59) 


or 


which  means  the  sliding  friction  moment  due  to  the  weight  W  is 

-*•  0)     ^  0) 

M._  =  y_WR_  [1  -  (-Ji)2]  _£  (61) 

sg     s   s        (jj  ^j^ 


The  torques  acting  on  the  obturator  are  the 
summation  of  the  torques  due  to  friction,  both  viscous  and 
sliding  (aerodynamic  or  gravitational) .   The  general  form 
of  the  torques  acting  to  oppose  the  motion  of  the  obturator 
in  terms  of  total  moment  M^  and  the  unit  vector  1   along  cu   is 

G  =  -M^l^  (62) 

which,  when  broken  into  the  various  contributors  becomes, 
for  the  projectile  in  flight. 


Gp  =  -(M^+Ms^)l3. 


(63) 
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TABLE  I 

COMMON  VALUES 

^o 

= 

87.25° 

*o 

= 

0 

K 

= 

0 

^o 

= 

0 

^o 

= 

0 

p 

= 

7669.6  Kg/m^ 

Rs 

= 

7.899  mm 

^p 

= 

4. 72  mm 

c 

z= 

3.125x10"^  N-m-s^ 

A 

= 

2.216x10"^  N-m-s^ 

X 

= 

0.41 

^f 

= 

1.917x10"^  Pa-s 

^s 

= 

0.35 

h 

=: 

5.08x10"^  m 

S 

= 

12.38x10"^  N-m-s 

AERODYNAMIC  MODEL              GRAVITATIONAL  MODEL 

^ 

= 

9.96  mm                W         =   0.07987  N 

\ 

= 

76.2  mm               E(eo,7r/2)  =  1.0063 

m 
P 

= 

0.141  Kg               cOp        =  200  sec"^ 

^P 

= 

8.58x10"^  N-m-s^ 

M 

1 

= 

3.0 

Y 

= 

1.4 

P 

1 

= 

0.099  MPa 

tii_ 

= 

12042  sec~^ 
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and  for  the  experimental  model, 

G  =  -CM  +M   )1  (64) 

m      V  sg  r 

These  may  easily  be  transformed  into  the  body- fixed  coor- 
dinate system  for  solution  of  the  equations  of  motion. 

For  solving  Eq.  (17) ,  a  numerical  integration  sub- 
routine [25]  was  used  and  computer  plots  of  6  vs.  time  were 
obtained.   A  copy  of  the  calling  program  is  included  in 
Appendix  A.   Figures  13  and  14  show  plotted  output  for  the 
input  parameters  shown  in  Table  1. 

C.   DEVELOPMENT  OF  AN  APPROXIMATE  SOLUTION 

-  Observations  of  the  nature  of  the  exact  solution  (Figs. 
13  and  14)  indicate  that  the  response  of  the  obturator  orien- 
tation (9)  to  the  applied  torques  (determined  from  Table  1 
parameters)  is  similar  to  that  of  a  linear  damped  second-order 
system  with  a  combined  step  and  ramp  input.   That  is,  the 
relationship  closely  follows  a  relationship  of  the  form 

^n  ^         d-c')^    '^n     "  ^ 

-  D  sin[a)n(l-C^)  t  +  (j)^^]  }  (65) 

in  which  the  parameters  are: 

D   =  step  height  (negative) 

R   =  ramp  slope  (negative) 

cOjj  =  natural  frequency 

5   =  damping  ratio 

^D'^R  ~  phase  angles 
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These  may  be  determined  by  means  of  a  solution  to  the 
linearized  equations  of  motion  under  the  assumption  of 
small  applied  torques  [  (M/Cco  ^)<<1].   This  has  been  accom- 
plished^  by  a  perturbation  of  the  steady  state  solution  given 
by  Eq.  (18)  for  the  case  of  no  applied  torques. 

The  first  step  in  the  approximation  is  to  rewrite  Eq.  (18) 
to  read 

X\l)   cose   +  (X+l)i   =  0  (66) 

^o    o         o 

where  ^p  ,  ^  and  9  are  constants.  The  dependent  variables 
of  the  general  form  of  the  governing  equations  (Eq.  17)  are 
written  as  sums  of  the  unforced  values  and  the  perturbation 
values.   Thus 


(67) 


The  general  form  for  the  applied  torques  on  the  obturator 
in  the  experimental  model  (gravity-driven  sliding  friction) 
is  given  by  Eq.  (64)  in  which  the  moment  term  is  composed  of 
a  fluid  (My.)  and  a  sliding  friction  (Mgg)  term.   When  both 
are  expanded  from  Eqs .  (42)  and  (61),  Eq.  (62)  becomes 

->.  OJ  ^        CO 

G  =   -    {Soj^  +    y3WR3     [l-(3^)=^]     }    ^  (68) 

r  r 


2 

The  theoretical  development  of  the  linear  approximation 
is  not  complete.   The  analysis  presented  here  therefore 
represents  a  preliminary  report. 

€9 


The  relative  velocity  co   from  Eqs .  (13,  35  and  36)  has  com- 
ponents 


oj    =  (il;-(i)  )sin9sind)  +  0cosd) 

ir  "^   p 

0),^  =  (ip-o)-)  sin9cos(j)  -  Bsintj) 
*        ^  (69) 

'^   y.  -  (i-oj  )cos9  +  4> 

3  i  ir 


To  the  first  approximation,  ^^-ui      (^q=«jo   for  the  case  of 
no  applied  torques) .   With  this  approximation,  and  the 
application  of  the  perturbed  variables,  Eq.  (67), 

u)   =  9cosd) 
ir   ~    o 

0)   =  -0sin(}) 
2r         o 

•  •         • 

0)   =  i|;cose^  +  (J)^  +  (J)  (70) 

3r   ~    ^  <->   ~ 

The  only  term  of  zeroeth  order  in  Eqs.  (70)  is  (^^   so 
that,  to  the  first  approximation, 

bi      -   ui        -  U^l  (71) 

The  modulus  is  taken  in  order  to  preserve  the  positive  sense 
of  0)^.   Under  the  same  approximations  as  above, 

'^x  =  *oSin<|)^sin^Q  (72) 

Substituting  Eqs.  (71  and  72)  into  Eq.  (60)  for  the 
sliding  friction  moment  arm,  Tg/Rg  becomes 


70 


-3.  =   Il-CsineQsinil;^)^]  (73) 


R 

s 


For  the  purpose  of  the  linear  approximation,  the  mean  value 
of  this  moment  arm  is  required.   Integrating  to  find  the 
average  moment  arm  gives 


r_  =   1   r   (r  /R  )  dt  (74) 

R     T     g  s 


where  T  is  the  period  for  one  revolution.   From  the  approxi- 
mation, however, 

T  =  ^  =  ^^^  (75) 

^o     ^o  I 

and  because  of  symmetry,  the  average  is  valid  over  one-quarter 

period,  therefore 

tt/4 
^  =  ^     f      (l-sin^9^sin^  ^  ) ""  dt  (76) 

o 

By   letting   u  =    27Tt/T   and   du  =    (27r/T)dt,    Eq.     (76)    becomes 


ir/2 


h 


I  =    -2     f        (l-sin^e^sin^u)    du  (77) 

o 


71 


Letting  EC9  )  denote  a  Complete  Elliptic  Integral  of  the 
Second  Kind  with  modular  angle  9  ,  Eq.  (77)  may  now  be 
written  as 


Solutions  to  Eq.  (78)  are  tabulated  for  various  values  of 

The  expression  for  the  applied  torque  vector  (Eq.  68) 
may  now  be  written  as 


G  =  -[SUI  +  U3WR3  I  EOq)].^  (79) 


With  these  approximations,  the  applied  torque  terms 
appearing  in  Eq.  (17)  are  given  by 


G         G 

1  2         -M 

^smct)  +  ^-coscj)  =  ^^(ip-00p)sin9-O 


G         G  * 

-lcos(})  -  -J-sin*  =  1^9  -  -I^ 


G  .  .  .        i 

_i  =  ZM   [(^_^  )co3e  +  <j,]  ^      -M    A  .  z^        o  (80) 

^     ^r       P  C\U       °    C  |i  I 


where  M  is  given  by 

M  =  S|i  I  +  y  WR   I  E(9^)  (81) 

O        S    S  7T      O 

and   a)p  is  the  projectile  angular  velocity. 

72 


Equations  (66) ,  (67)  and  (80)  are  combined  with  Eq.  (17) 
After  simplification,  and  retaining  only  the  first  order 
(linear)  terms,  the  results  are: 


i  =  -(ii^Qcote^)  e 

..  •  'Mo 

J  =    i^l)^cscQ^)Q   -  ^  — — 


(82) 
(83) 


-M     *  *  ,  •  *  •  -, 

e   =       ,  .       9    -    il^^sineo    [(j)(l+X)-eXi|^^sine^   +    ii;XcosG^]  (84) 

aI^qI"-  ~  ~°         o~  o 


These  equations  are  then  integrated  under  the  following 
initial  conditions 


^ '      ^o    p 


({)(0)  =  0 


6(0)  =  e 
9(0)  =  0 


} 


\l)   =   0 


*^°^  =  -*o  =  ITT^o^°^^ 


9(0)  =  0 
9(0)  =  0 


Eq.  (82)  yields 


i  =  -i^i)^cote^)e 


(85) 


From  Eq.  (8  3) 


M. 


(j)  =  -4)   +  (ij;  cos 9^)9  -  iit-x 

~   C  1 4) 


*, 


(86j 


By  substituting  the  expressions  obtained  in  Eqs.  (85)  and  (86) 
into  Eq.  (84) ,  the  differential  equation  for  9  becomes 


.73 


e  =  ^i^  -  i  sine^Ei^e  -  (i+x)  (i^  +  ^  t)  ]  (87) 


or,    in   a  more    familiar   form. 


-»-  +  ^=-^  e  +  e  =  D  + 


2     .  .    .     .         ^    .    Rt  (88) 


whe 


^®   %  ^   ^o  ^   %'    ^0/1*^0!    =   "^   ^^ 


C=  ^ 


2Ai(j)o|i|io 


D  =  il±Al4,^sineQ  =  -  ^inae^ 


R  =  -  lilAL  M3ine    =  -  -^ 


^ 


CO 
o  <t>oC 


Expansion  of  the  expressions  above,  using  Eq.  (81)  shows 
that  the  damping  coefficient  is  given  by: 


or 


C  =  }    ,.   [S|i  I  +  y  WR  |e(0o)] 

2aUoUo  s   stt   o 


C  =  —J-  [S  +  ygRgW  i  _4A±1^ E(eo)J 

2AiiJo  ^  x^o^oseo 


and  the  ramp  slope  (rate  of  decay)  becomes 


7.4 


X 
R  =  ^-  TA+iy'^o^'-^2®o)^  (90) 


The  integral  of  Eq.  (87)  with  the  initial  conditions 
specified  [0(0)  =  6(0)  =  0]  is  Eq.  (65)  where 

h 
*R  "  ^^D^   2tan-^  JllCil  ^g^j 


For  large  values  of  C(i)>.t/  the  ramp  decay  dominates  the 
response  and  the  solution  may  be  further  approximated  by 

e=e+e  =  9+D  +  ~(a)  t-2c) 

o    ^    o       coj^   n 


The  time  to  nutate  through  a  known  angle  9  may  be  found 
from 

e-e  -D   ,    ,    (9  -9) 
t  =     Q    +  2C  ^  2£  _  _o _  D_  (92) 

R     a)„    aj„      R      R 
n    n 

where  r,,    R,  D  and  oo^^  are  prescribed  in  Eqs .  (88) -(90). 

Figure  15  is  a  plot  of  6  vs.  time  comparing  the  results 
of  the  exact  solution  from  Fig.  14  with  those  of  the 
approximation  developed  here.   The  physical  constants  and 
initial  conditions  are  identical  for  the  two  cases. 

D.   THE  EFFECT  OF  GEOMETRY  ON  OBTURATOR  RESPONSE 

In  the  discussion  to  this  point  a  primary  consideration 
has  been  to  model  the  external  forces  acting  on  the  obturator, 
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The  effect  of  geometry  has  not  been  discussed,  but  is  a  most 
significant  contributor  to  the  response  of  the  obturator. 

C  and  A  are  functions  of  mass  density  and  physical  di- 
mensions.  By  varying  the  material  composition  (density) 
the  response  will  vary.   Utilizing  the  physical  dimensions 
of  the  obturator  in  Table  i,  the  obturator  volume  may  be 
calculated  by  subtracting  the  volumes  of  a  right  circular 
cylinder  of  radius  r  ,  length  2RgCOs6s  ^^<^   the  two  spherical 
caps  at  the  ends  of  the  hole  from  the  volume  of  a  sphere  [26] 
Thus  the  volume  of  the  obturator  is  given  by 

A  O 

V  =  tttR„^-t^(R  -R^cose)  ^  (2R  +  R^cosB)  -  27rr„^R^cos3 

■joJSS  ss  ps 

Collecting  terms  gives 

2 


V  =  2TrRsCOs3[Rs'(l  "  ^^f-^)  "  r^'] 


By  knowing  the  mass  density,  the  weight  of  the  obturator 
[for  use  in  the  torque  expression,  Eq,  (61) ]  may  be  obtained, 
Therefore  as  the  density  varies,  the  values  of  C,  A  and  W 
vary.   Fig.  16  is  a  family  of  curves  of  6  vs.  t  for  varying 
density  with  all  other  physical  properties  being  those  of 
Table  1. 

There  is  however  a  limiting  case  for  the  size  of  the 
hole.   If  the  hole  becomes  too  large  (assuming  same  size 
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projectile  hole) ,  there  may  not  be  sufficient  obturator 
surface  area  to  block,  the  passage.   The  limiting  relation- 
ship is 


"s      ^ 


This  limit  is  purely  theoretical  and  is  not  a  feasible  design 
parameter.   Engineering  judgement  in  the  selection  of  realis- 
tic values  is  required  so  that  the  projectile  will  meet  the 
ballistics  criteria  and  retain  its  integrity  while  undergoing 
the  forces  associated  with  the  entire  firing  evolution  from 
chambering  to  target  impact. 

The  effect  on  the  obturator  response  due  to  increasing  the 
hole  diameter  is  shown  in  Fig.  17  .   The  inertial  contribution 
to  the  response  (coupled  term)  is  increased  by  increasing  C/A 
and  therefore  X;  however  the  friction  torque  (Mg_)  contribution 
is  lessened  by  the  reduction  in  the  mass  of  the  ball. 
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III.  EXPERIMENTATION 

A.   EXPERIMENTAL  APPARATUS 

An  experimental  system  was  designed  to  simulate  the 
spinning  BOT.   The  apparatus  consisted  of  a  compressed-air- 
driven  spin-up  rig,  optical    timing  mechanism,  air  manifold 
and  associated  piping,  tubing  and  electronics. 

The  obturator  was  fabricated  from  a  standard  5/8- in. 
diameter  chrome -steel  (52100)  bearing  ball.   The  ball  was 
annealed  to  allow  machining  and  then  bored  along  a  diametrical 
axis  (Fig.  18) .   The  bored  ball  was  then  mounted  in  a  bakelite 
metalographic  specimen  mount.   The  mounted  ball  was  placed  in  a 
milling  machine  and  a  flat  was  machined  in  the  bakelite  at  a 
specified  angle  relative  to  the  axis  through  the  hole  in  the 
ball  (Figs.  19  and  2  0) .   This  flat  was  then  used  as  the 
polishing  plane  for  metalographic  specimen  preparation.   A 
small  flat  spot  was  polished  on  the  obturator  to  provide  a 
highly  reflective  surface  at  a  known  orientation  relative  to 
to  the  z-axis  of  the  obturator  (Fig.  21) . 

After  removal  from  the  bakelite  mount,   the  obturator 
was  placed  in  a  three-piece,  lucite  housing  (Fig.  22).   The 
mating  ends  of  the  two  hollow  inner  cylinders  were  each 
machined  with  a  5/8- inch  end  mill  to  a  depth  of  approximately 
5/16-inch.   When  mated,  a  spherical  cavity  was  formed  to 
accomodate  the  obturator.   The  third  cylinder  was  press  fit 
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Figure    18. 


Boring  of  hole    along   a  diametrical   axis    of 
the   annealed  bearing  ball. 


Figure  19 . 


Machining  of  flat  in  metalographic  specimen 
mount. 


82 


,  -s*^  s^#*i«^*^*s  ".i^'"*^^'^, "  »'V~>*'f*^''--T^ 


Figure  20.   Machined  specimen  mounts  with  varying 
angles  (6  ) . 


Figure  21.   Polished  spot  orientation  on  obturator  and 
relationship  to  9  . 
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Figure  22. 


Obturator  and  components  of  lucite  obturator 
housing  with  mating  ends  of  inner  cylinders 
machined  to  house  the  obturator. 
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Figure  23.   Assembled  obturator  housing  with  aluminum 
end  pieces. 
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over  the  others  after  the  obturator  was  inserted  to  insure 
alignment  of  the  inner  pieces  and  provide  rigidity  (Fig.  23). 
The  ends  of  this  lucite  assembly  were  then  press  fit  into 
aluminum  end-pieces  similar  to  those  in  Fig.  23 .   These  end- 
pieces  served  as  the  shaft  for  the  bearings  and  one  also  served 
as  the  prime  mover  (bucket  wheel)  for  the  apparatus.   The  shaft 
rotated  in  two  ball  bearings  mounted  in  aluminum  pillow  blocks 
aligned  on  a  rigid  pedestal.   The  prime  mover  was  a  bucket 
wheel  machined  from  a  solid  aluminim  disk  (^ig.  24) , 

The  bucket  wheel  was  driven  by  compressed  air  supplied 
from  an  installed  system  through  an  air  filter  to  a  manifold 
and  then  through  one  of  two  Model-10  Kendall  pressure  regu- 
lators.  From  the  regulator,  the  air  passed  through  a  flexible 
tube  to  a  tee,  each  leg  of  which  supplied  a  nozzle.   These 
nozzles  were  mounted  opposite  one  another  on  the  pillow  block 
in  such  a  way  as  to  allow  the  air  jet  to  impinge  upon  the  bucket 
wheel  to  cause  rotation  (Fig.  25) .   The  speed  of  rotation  was 
sensed  by  a  Bentley  Nevada  Proximitor,  Model  3100N,  which  was 
mounted  above  the  bucket  wheel  to  detect  the  passage  of  each 
point  on  the  wheel.   The  sensor  was  supplied  from  a  LAMBDA 
Regulator  Power  Supply  Model  LP413FM  by  18VDC.   The  pulses 
generated  by  the  sensor  were  counted,  averaged  and  displayed 
as  a  frequency  by  a  Monsanto  Programmable  Counter-Timer  Model 
HOB. 

While  the  spin-up  rig  was  being  brought  up  to  the  desired 
speed,  an  air  jet  from  a  nozzle  mounted  rigidly  on  the  pedestal 
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at  the  end  opposite  the  bucket  wheel  held  the  obturator 
fixed  to  the  spinning  rig.   This  jet  was  supplied  from  the 
common  air  manifold  and  the  other  pressure  regulator  also 
through  flexible  tubing.   The  air  jet  passed  through  a  hole 
in  the  plunger  of  a  Rocker  Solenoid,  R.S.  No.  10-207.   This 
hole  was  aligned  with  the  hole  through  the  "projectile"  to 
allow  the  jet  to  impinge  upon  the  obturator  (Fig.  26) .   When 
the  switch  at  the  far  right  of  Fig.  26  was  activated,  standard 
115VAC  was  applied  to  the  Rocker  solenoid  and  removed  from  a 
normally  closed  ASCO  Solenoid  Valve  in  the  air  supply  line  to 
the  jet.   The  plunger  retracted  from  the  position  shown  in 
Fig.  26  and  the  solenoid  valve  closed.   The  retraction  of  the 
plunger  performed  three  functions.   It  first  caused  a  pulse 
to  be  generated  by  another  Proximitor  mounted  next  to  the 
nozzle  and  powered  from  the  same  18VDC  supply.   This  pulse 
started  the  timer  function  of  a  second  Model  HOB  Counter-Timer, 
The  plunger  also  covered  the  nozzle  outlet  by  misaligning  its 
hole  with  that  of  the  nozzle.   This  removed  the  restraining  in- 
fluence of  the  air  jet  on  the  obturator  and  prevented  any  air 
remaining  in  the  supply  line  from  impinging  on  the  obturator. 
And  lastly,  in  the  fully  retracted  position,  the  machined  and 
polished  end  of  the  plunger  was  positioned  opposite  the  pro- 
jectile hole  (Fig.  27) .   When  the  released  obturator  nutated 
through  enough  of  an  angle,  the  beam  of  a  SPECTRA  PHYSICS  MODEL 
132  LASER,  MODEL  NO.  3187,  passed  through  the  projectile  and 
was  reflected  by  the  polished  end  of  the  plunger  into  a  light 
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sensitive  diode  (Fig.  2  8) .   The  pulse  emitted  by  the  diode 
triggered  the  stop  channel  en  the  timer  and  the  elapsed  time 
was  displayed.   Figure  29  shows  the  entire  experimental  set-up 
with  air  flasks  in  the  background. 

B.   EXPERIMENTAL  PROCEDURE 

The  obturator  was  positioned  in  the  spin-up  rig  in  such 
a  way  that  the  laser  beam  was  reflected  off  the  polished 
spot  back  to  the  target  mounted  on  the  laser  (Figs.  29  and 
30) .   The  center  of  the  target  is  the  location  of  the  beam 
and  the  obturator  was  adjusted  to  place  the  reflected  spot 
on  the  center  of  the  target.   Thus  the  orientation  of  the 
z-axis  through  the  obturator  was  known  relative  to  the  Z-axis 
(laser  beam)  through  the  projectile.   This  is  the  initial 
value  9q. 

The  solenoid  plunger  was  checked  in  position  with  the 
hole  aligned  with  the  nozzle  (Fig.  26)  and  then  air  was 
admitted  to  hold  the  obturator  in  place.   The  pressure  was 
then  increased  to  the  turbine  to  set  the  apparatus  in  motion. 
Figure  31  illustrates  the  circle  traced  by  the  reflected  spot 
(Fig.  30)  on  the  target.   This  circle  facilitates  alignment 
of  the  laser  to  insure  the  beam  was  centered  through  the  pro- 
jectile.  (It  can  be  seen  in  Fig.  31  that  the  beam  was  slightly 
off-center.)   Alignment  was  accomplished  by  positioning  the 
beam  from  side  to  side  for  horizontal  alignment,  then  raising 
or  lowering  the  spin-up  rig's  mounting  platform  with  adjusting 
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Figure  30.   Laser  reflection  on  the  target  from  the 
polished  spot  on  the  obturator. 


Figure  31. 


Circle  traced  by  the  reflected  laser  while 
the  apparatus  is  spinning. 
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screws  for  vertical  alignment.  When  aligned,  the  circular 
trace  of  the  reflection  exactly  follows  the  circles  of  the 
target. 

Once  the  alignment  had  been  checked  and  the  obturator 
positioned  at  the  known  9  ,  the  pressure  was  increased  to 
the  turbine  nozzles  until  the  desired  spin  rate  was  achieved. 
The  regulator  was  then  adjusted  so  as  to  maintain  the  spin 
rate  at  ±5  Hz  on  the  digital  display  or  approximately  ±1.75 
rad/sec.   (The  actual  tubular  projectile  spin  rate  is 
approximately  1.2x10**  rad/sec).   Once  the  desired  rate  was 
achieved  and  noted,  and  the  timer  checked  and  reset  if  neces- 
sary, the  switch  (Fig.  26)  was  activated.   This  secured  the 
air  to  the  nozzle  holding  the  obturator  and  retracted  the 
solenoid  plunger  (Fig.  2  7) .   The  obturator  was  released  and 
began  to  move  relative  to  the  spin-up  rig.   When  the  angle 
9  reached  the  value  at  which  the  laser  beam  could  pass  through 
the  obturator  hole  iQ=Qf^-^=B>)  ,  the  beam  struck  the  polished 
plunger  and  was  reflected  into  the  light  sensitive  diode 
housing  (Fig.  28) .   The  reflected  beam  striking  the  diode 
caused  the  diode  to  emit  a  pulse  which  turned  off  the  timer 
gate.   The  elapsed  time  from  plunger  retraction  to  release 

the  obturator,  to  the  obturator  nutating  to  0   was  dis- 

^     cr 

played  on  the  timer  display  to  the  0.0001  second.   This  time, 

t   ,  was  recorded. 
cr 
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C.   DATA  REDUCTION 

The  only  data  reduction  required  was  the  conversion  from 
Hz  to  Radians  per  second.   The  frequency  read-out  on  the 
digital  display  indicated  the  number  of  teeth  on  the  bucket- 
wheel  passing  the  sensor  in  one  second.   Since  there  were 
eighteen  teeth  on  the  bucket  wheel,  the  conversion  was 
simply 


f  f    TT 

0)^  =  -^2  X  277  =   ^Q  sec"^ 
""P     18  9 


Appendix  B  is  a  listing  of  the  raw  data  and  the 


corresponding  values  of  oo  , 


9^6 


IV.  CORRELATION  OF  RESULTS 

Figures  32  through  3  7  show  0  vs.  time  for  various  values 
of  (o_  from  40  to  1000  sec""^.   At  high  and  low  values  of  co,-,, 
9  approaches  zero  more  slowly  than  at  the  intermediate  values. 
Figure  38  is  a  plot  of  Eq.  (92)  and  the  experimental  data. 
The  plot  of  Eq.  (92)  uses  the  physical  qualities  from  Table  1 
and  y  =0.35  over  the  range  100-^-^100  0  sec"^.   It  is  clear 

S  tr 

that  the  experimental  data  falls  below  the  approximation  curve 

at  high  o)   and  rises  sharply  above  the  curve  at  low  oo  . 

In  order  to  determine  what  parameters  may  effect  the 

position  of  the  approximate  curve,  sensitivity  checks  were 

run  on  the  quantities  in  which  the  most  uncertainty  existed, 

namely  Ps  ^^^  ^o  (®"^®^  though  Qq   could  only  vary  by  ±1°)  . 

Figures  39  and  40  indicate  the  sensitivity  of  t   to  variations 

in  u   and  6_  respectively.   Small  variations  in   0   are  not 
s      o    '^       -^  o 

discernable  within  the  range  of  scatter  of  the  data  (Fig.  41). 
Variations  in  Pg  have  a  much  more  pronounced  effect.   From 
Fig.  42  it  may  be  seen  that  the  data  at  higher  values  of  co 
correlate  well  for  a  y  =0.55. 

It  has  been  noted  that  the  approximate  solution  is  based 
upon  the  assumption  M/Coo  ^<<1.   At  values  of  oo  (<200  sec~M/ 
this  assumption  is  no  longer  valid  and  therefore  calculations 
and  correlation  of  the  data  to  the  approximate  solution  will 
not  be  conducted  for  a)„  200  sec"  ^ .   Also  at  low  values  of  co_, 
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t   appears  to  approach  <»  almost  assymptotically .   Because 
9o-90°  (87.25°),  the  obturator  may  be  near  a  metastable  state 
at  low  0)   such  that  the  coupled  terms  in  the  equations  of 
motion  have  little  effect.   The  amplitude  of  oscillation  at 
low  ojp  is  also  relatively  small. 

In  general,  for  large  values  of  ^'^j^t  (Eq.  65)  and 
neglecting  the  effect  of  viscosity,  t  becomes 


t  =  A±l  M +   ^"^P   [(9.-6)-  ^in9^]  (93) 

^   Aoip^cose^    Msin9^     o      2     o 


where  M  is  the  gravitational  contribution  from  Eq.  (81). 
This  may  be  written  as 

T 

t  =    I   +  T  03^  (94) 

3       2  P 

with  T   and  T  being  constants  determined  from  the  known 

1         2 

physical  parameters  (Table  1) .   Utilizing  values  of 

Oq  =  87. 25°,  9  =  e^^  =  sin"^  (r/Rg)  =  36.734°  and  the  measured 

weight  of  the  obturator,  V7  =  0.079  8  N, 

M  is  found  to  be 

M  =  6.34x10"**  y  N-m 

Therefore, 

T   =  2.039x10^  jl 

and 

T   =  3.03x10"'*  /  y  (95) 

2 


'109 


thus 


t    2.039x10^   -  , 


^   ~   3.0  3x10"*   ^ 


p 


But  for  large  values  of  co  , 


3. 03x10"'*  03 
t  =  £  (96) 


A  least  squares  fit  of  the  data  Cc0p^35.97)  to  Eq.  (96)/ 


such  that  t  =  Toop,  gives 


St. 
T     =        ^ 


2 
2 


2         ^tia)pi 

Substituting   in   the   data  points    (denoted  by    '*'    in  Appendix  B) , 

Ztj_^      =      2.1046884    sec^ 

Eojpi^    =      2.8413714xl0^sec"^ 

Zti^O)    .  ^    =    2.430257x10^ 
^      pi 

therefore 

T      =    8.6603508x10""*    sec^ 
Since  the   standard  deviation    for  T   is    given   by 


=  ^=-      zt..^ 


a^'  o^ll\  1 


110 


where 


0^  =  ^  C^^pi'--r2^tia)p.) 

the  following  results  are  obtained: 

Q^   =  1.4659547x10^ 

a,p2    =   3.9181053x10"^° 

a„     =    1.97742x10"= 
T 

-i     =    0.023 
T 

The  data  therefore,  deviates  from  the  least  squares  value 
by  only  2.3  percent. 

Using  the  value  calculated  for  T   in  Eq.  (95) , 

2 

_    3.03x10"'^ 

U   =  rr  =0.35 

8.66x10  ^ 

with  EC87.25°)  =  1.00458, 


£_  =   2  (1.00448)  =  0.63954 


thus 

0.35 
Ps  =  0.63954  =  0.5473 


111 


with  this  value  of  y  ,  values  of  t   were  obtained  from 

s  cr 


._  1 


the  exact  solution  over  the  range  12  5<.(A)p<.70  0  sec"   and 
plotted  with  the  approximation  of  t^^  from  Eq.  94  and  the  data 
points  (Fig.  43) . 
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V.  CONCLUSIONS 

The  response  of  the  obturator  (9)  as  a  function  of  time 
to  applied  torques  was  found  to  be  that  of  a  rigid  body  of 
revolution  about  its  mass  center.   Euler's  Modified  Equations 
of  Motion  in  terms  of  Euler  Angles  describe  the  response  to 
these  applied  torques.   The  response  is  predictable  and  was 
compared  for  various  combinations  of  physical  parameters.   By 
varying  density,  hole  size  and  angular  velocity,  various  re- 
sponse  times  may  be  obtained.   It  was  found  that  in  each 
case,  a  minimum  response  time  was  attainable  for  a  given  value 
of  a)p.   For  the  case  of  the  experimental  apparatus,  this  mini- 
mum occured  at  oa   =  150  sec~^. 

P 

The  analytical  model  and  a  linear  approximation  to  this 
model  were  used  with  the  parameters  of  the  experimental 
apparatus  and  close  agreement  to  the  data  was  achieved  at 
w^>200  sec"^.   Below  ojp  =  200  sec~^  the  exact  solution  appeared 
to  be  higher  and  the  approximate  solution  lower  than  the  experi- 
mentally observed  values  of  t^r*   This  is  explainable  in  the 
approximation  since  it  was  based  upon  the  assumption  that 
M/a)p^C<<l.   At  low  values  of  Wp,  this  assumption  does  not  hold. 

The  trend  in  the  exact  solution  at  low  ojp  is  to  approach 
^cr  ^  °°   almost  assymptotically.   This  trend  is  observable  in 
the  experimental  data  as  well  and  may  be  explained  by  the 
metastable  initial  position  of  the  obturator  9^-90°  (9o=90°  is 
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a  stable  position) .   The  effect  of  the  coupled  terms  in  the 
equations  of  motion  is  small  at  low  w  ,    therefore  the  initial 
effects  are  not  as  great  and  the  amplitude  of  oscillation  of 
0  is  insufficient  to  cause  departure  from  this  metastable 
position.   Since  the  actual  projectile  operates  at  a)p-12,000 
sec~^,  high  values  of  a)p  are  of  interest.   The  theory  agrees 
well  with  the  data  above  a)p=200  sec*"^  and  therefore  the 
correlation  to  the  data  was  made  at  these  higher  values. 

Because  the  linear  approximation  agreed  closely  at  higher 
ujp,  a  least  squares  fit  of  that  approximate  equation  to  the 
experimental  data  (a)p— 235.97  sec~M  lead  to  the  calculation 
of  a  higher  sliding  friction  coefficient  than  had  been  used 
previously.   The  assumed  value  was  yg=0.35  for  laminated 
plastic  on  steel  [27].   This  was  found  to  be  too  low  by  the 
least  squares  fit.   A  sensitivity  check  of  the  approximate 
solution  for  t^j-  to  variations  in  Pg  also  indicated  that  0.35 
was  too  low  (Fig.  39).   Therefore,  a  method  for  determining 
the  sliding  friction  coefficient  from  the  experimental  data 
and  the  approximate  solution  was  discovered. 

By  knowing  initial  orientation  (9q) ,    projectile  spin  rate 
(ojp)  and  the  physical  dimensions  and  material  properties  of 
the  projectile,  the  response  may  be  predicted  in  a  gravity 
environment  with  reasonable  accuracy. 

The  aerodynamic  model  appears  to  give  reasonable  results 
within  the  scope  of  this  study.   The  complex  forces  acting  on 
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the  obturator  in  flight  are  extremely  difficult  to  model 
theoretically  and  therefore  they  were  simplified  greatly. 
The  confidence  gained  in  the  gravity  model,  however,  leads 
to  the  conclusion  that  once  the  complex  aerodynamic  forces 
are  better  understood,  the  motion  of  the  ball  within  the 
flight  model  will  be  predictable. 
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VI.  RECOMMENDATIONS 

Further  experimentation  should  be  undertaken  with 
variation  of  other  parameters  (hole  size,  obturator  material) 
to  obtain  more  confidence  in  the  analytical  model.   Wind 
tunnel  tests  of  the  BOT  should  be  conducted  in  order  to  more 
closely  approximate  the  actual  forces  on  the  obturator  in 
flight,  and  to  determine  the  position  (6)  at  which  the  de- 
tached bow  shock  is  swallowed.   This  would  greatly  enhance 
the  predictive  capabilities  of  the  model  and  thus  allow 
coupling  to  a  trajectory  model  for  predicting  flight  char- 
acteristics.  This,  in  turn,  could  be  coupled  to  a  numerical 
optimization  routine  for  optimizing  the  BOT  design  for  the 
desired  performance  criteria. 
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APPENDIX   A   -    COMPUTER   CODE 

C  «    BALL    QBTUPAT'-D    TUBULAR    PRQJ>:CTILn    • 

C  .              »  BALL    MOTION    ANALYSIS                   * 

C  ♦  ♦ 

C  *  LT   J,    W,    BLCDMER    lit    USN                ♦ 

C  *  * 

r  *****««*««**♦****  *«»*^*:»t  ************* 

OIMtfNSICN    XOO)  ,XnnT(30»,C(  15) 
C  C(l)=    (SIN(R[:T/ )  1**2,    Wh-R""    Bl^TA    IS    GIV^N    BY   ARSIN(P.P /RO) 

C  C(2»=    L/iMBCA    =     (C-A)/A 

C  C(3)=    «-kCOYN*'MIC    np     GRAVITATIONAL    «;.LI'DINIG    FRICTICN   CCNST4NT. 

C  C(4)=    TCTAL    Kir.r-TTc    ^Hi-.RC-y    Of    CCMBIN'^D    PROJlCTIL"    AMD    OBTURATOR. 

C  C(5>=    MIMMLH    yc»"-N-^    '1=    IN^^cTIA,    A. 

C  C(6)=    .--lAXLVLlM    f-CMrNT    OF     I^.«=RTIA,    C. 

C  C(7)a    MC^'£NT    CF    IM^oyxA    C^    A    J-OLLOW    CIRCULAR    CYL  INO!;^ 

C  APPROXI^^TIMG    THAT    ZF    THt    PROJCCTI  Lf^, 

C  C(8)=    FLUID    FSICTICN    CONSTANT,    S. 

C(10)  =  U 

c  • 

c 

C»*****»CALL  INTEG2  SIVULTANCCUS  OIFF.  EON.  SOLVER 
C 

c 

I  CALL  INTcG2<T,X,X00T,C) 
C 

c 

C*****«*OtFINr  ELLEF  ANGLTS  Th'TA,  PHI  AND  PSI  AND  1ST  OERIVATIVTS 

I 

TH=X(1) 

X0CT(1)=X(2) 

PS=X<3) 

X0CT(3)=X{A) 

PHsX(5) 

XCCT(5)=X(6) 

i    . 

C*******D£:FINf'  TRIGCNCMSTRIC  FUNCTIONS  OF  EULFR  ANGLES 
C 

c 

STI-  =  SIN(TH» 
CTh  =  COS(Tl-) 
SPH=SINiFH) 
CP»-  =  COS(P»- ) 
SP«=SIN(PS) 
CP£=CJSIPS) 

c 

C*******Cf:FINu    TI-ETi-DOT,    PHI-CCT    AND    PSI-DOT 
C 

TI-C=XDGTC1) 

PSD=X0rT(3) 
PH0=XDOT(5) 

c 

C****^***CALCULAT3    ANGULAR    V^LGCITIcS    WITH    R'SPSCT    TO    BODY- 

C  FIX£D    AXcS    <^,    B,    ANO    C    Ri  SP~CTI  V"LV. 

C 

OMl  =  PSD«STH*SFH  +  THn*-rpH 

aW2=PSD»ST>-*CFH-TH0*;$PH 

CN3  =  PS0*CTh*F»-C 
C 

C*******CALCULATi    ANGULAR    V?!LOC!TI!FS    ABOUT    INERTIAL    AXOS    X,    Y,    AND 
C  2    RcSPtCTIVflLY. 

C 

0^'X=GM1*  (CPH*rP3-SPH'«C'H*£PS)+C"2«(-SPH*CPS-CPH*CTH^SP3) 

*  +  CK3*(STK'-'?P5  ) 

0KY  =  3M1*(CPH*SP<"*SPH«CTH*CPS)  +  CM2*(-SPH»'SPS  +  CPH»CTH*CPS) 

♦  ♦0K3"^(-STI-*CPS) 
Q»'Z=CMl*STh*SPH4-Ci«2*STH-»CPH-fCM3*CTH 

C 

C**»***»*CALCUL4TF  TF""  ANGULAR  V"LCCITY  OF  TH:r  PROJECTILE 

C         FROM  CCNSiFVATi:^N  OF  r.N.RGY  PRIMCIPLcS. 

C 
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SPIN=0.0 

SPlN=SOrlT(  (2.*C(4>-C(5)*{CMl**2+aM2**2»-C(6)*0M2*«2)/CJ7») 
C 

C«-«**i.*»***CALCULATc    ANGULAR    VLOCITY   OF    BALL    REUTIVE    TO    TH?    PROJ^CTIL 
C  WITH    ?f:SP-:CT   TO    80CY-FIXED   AXTIS. 

c 

OMlRsaMl-SPlN*SPH«<;TH 
Qf2?l=CI»2-Sf  IN*CPI-*STH 
0HiR»C^«3-SFIN»CTH 

c 

C*»»»*«#*CALCULATr    fiiLATIV"    ANGULAR   V-LOCITI-S    WITH   RESPECT    TO 
C  THi    INIrtTIAL    PrF!  ^CNCr    PRAM^. 

c 

CHXR=af«lP=*(CFh*CPS-SPH*CTH*SPS»*0M2R*(-SPH*CPS-CPH-CTH*SoS) 
•♦Cf3R«(STH«SPS» 

QMYR=OMlR*(CPh*5PS+SPH*CTW*CPS)+OM2R*(-SPH*SPS+CPH*CTH«C«>S» 
♦♦CF3R*(-STH«CPS» 

OMZK=OMlR»STH*S?H+nM2R*STH*CPH+OM3R*CTH 

C***«'^***»CALCULAT?    fAGNITUO"    CF    RCLATIVT    Af^GULAR   veLCCITY    VGCTO^. 
C 

0^P=SQRT(CMF*C«lR+0M2R«Cf'2R-»CM3R*0M3R) 
C 

C********CALCULATE    TH=    MOMENT    fiRM    AT    WHICH    THE    FRICTIONAL    FORCES    ACT, 
C 

C  ARH=    MOMENT    ARM    AND   C(l)     TS    AS    0<;;FIN=D    SXCdPT    IN    THG    CASS    OF 

C  GRAVITY    WH«r:    C(I)     =    1.>. 

C  FRIC=   CGMBINcC    FLUID    AND    SLIDING    FRICTION    CONSTANTS. 

C 

ARF=S0RT(l.-C(l)*(ayXR/CKR|«*2-<l.-C(l) J*CQMZR/CMR)**2) 

FRIC=C(3)^AR^' 

AMCM=C(&  )*CI»'P+FRIC 

C****A*****CALCULATn  MC»«:NTS  ACTING  ON  THr  BALL  ABOUT  TH=  BODY-FIXED 

C  AXbS. 

C 

Gl'-AMQM  ♦uf' 1  ? /CO'R 
^      G2=-AMCK*C>2P/CVR 

C»**»**»**»CALCULA7^      ■LGM'^NTS    OF    MATRIX    OF    COEFFICIENTS »    Yl,    Y2    ANC    Y3 . 
C 

Y1=G1/C(  5)-PSD*f'^H0P':'^H*?PF-»>PI-C'«STH*CPH)+THD*PHC*SPH-C(2)* 
*(PSD*STh=»CPH-ThD*"5i^M)*(  P  5n*-C  TH  +  PHO  ) 

Y2  =  G2/C(5)-FSC*(~H-«'rTH*CPH-o»-D*STH*SPH)+PI-C*THC*CPH+C(  2)' 
*(PhO+PSO*C  ih  »»(P«:0  •SrH»':PH  +  THO-«CPH> 
Y3=G3/C(6) ♦PSC*">0*5TH 

c 

C***«*****»*CALCULAT.-:    ANGULAP    ACC!!L  =  R.ATI ON   T=iRMS    FOR    INTEGRATION,    2N0 

C  Of PIVATIV'S    OF    "ULTR    ANGL'S    WRT    TIMS. 

C 

XDOT(2)  =  CPh'Yl-S?H'=Y2 

XCCT(6)  =  -SFI-*CTH/STH*Yl-C^h*CPH/STH*Y2*Y3 

X0CT(4)  =  SPh/STH«Yl  +  CPH/S"fH*Y2 
C 

c 

C****^**CALCULATr    VALU'iS    FOR    THG    APPROXIMAT!:    SOLUTION 
C 

c 

THCRIT=.tAll 

IF(T.N'i.C.CCCCOCO)    GO    TO    25 
THNCT=X(1) 
•     0MP=X(4) 
FRTC=C<6) 
A=C(5) 
AMCa=C(2) 
A^CAP=AMCA♦1.     . 
E  N  =C  (  3  ^^ .  6  3  <:  f  A 
D  =  -AMnA/2.-»SIN(2.*THNOT) 
R=-AyCAP*:N=«SIN(THM3T)/  (     P'^C^O^'P) 
TCRI T=-l ./R*{ THNOT-THC  P ! ~ )-D/R 
WRIT'^(6,2M    TCRTT 
•       20    FCPMAr(5X,''fCP!~    =    •,"16.R) 

Zb7A=AMDAP/(2.«A*Af^nA*CNF*CMP*CGS(THNOT)  J*'M 
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IF(ZETA,Gr. 1.0)00   to    25 
DZtTA=SCKT(l  .-Z""TA*Z':ta) 

PHAS-=.^TAN{CZ':Ta/(Z-TA*(l.-2./(A.«ZnTA*Z*-TA-AMDAP))  )  » 
THFRaX=Th:NlOT*C+?/aMP^<-2.*?tTA)+l./ZKTA*SIN(PHASt) 
GO   TC    26 

25  THP.-<0X=Th:>>inT*D*R/GMP-(C>'P*T-2.'«Z[:TAI-H./  (EXP  (  ZtTA*OMP*T  )*OZ::TA  I* 
*SIN(OMP*T*CZu:TA*PHASi:) 

26  X(7)=!HPRCX 
X(5l=0MZ 
X(13)=SPIN 
GO    TO    1 
END 
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APPENDIX   B    - 

-    RAW    DATA 

READOUT 

TIME 

''p 
[sec-^J 

READOUT 

TII4E 

'"p 
[sec-M 

[HzJ 

Isecj 

[HzJ 

[secj 

776 

.274 

270.88* 

880 

.2717 

307.2* 

676 

.218 

235.97* 

920 

.3026 

321.14* 

679 

.2181 

237.02* 

920 

.2847 

321.14* 

573 

.2157 

200.01 

9  36 

.2816 

326.7* 

570 

.2079 

198.97 

983 

.2699 

343.1* 

574 

.2009 

200.36 

1001 

.306 

349.4* 

570 

.197 

198.97 

994 

.2972 

346.97* 

571 

.2031 

199.32 

1055 

.296 

368.26* 

475 

.2099 

165.8 

1063 

.2530 

371.06* 

479 

.1833 

16  7.2 

1151 

.315 

401.77* 

570 

.2179 

198.97 

1022 

.369 

356.7    * 

475 

.215 

165.8 

1175 

.378 

410.15* 

479 

.223 

167.2 

1258 

.3348 

439.12* 

470 

.2  30 

164.1 

1278 

.3968 

446.1   * 

428 

.1831 

149.4 

1255 

.4268 

438.08* 

425 

.1786 

148.35 

1221 

.3088 

426.2  * 

430 

.1659 

150.1 

1170 

.284 

408.4  * 

367 

.1896 

128.1 

1158 

.281 

404.2* 

363 

.2037 

126.7 

370 

.1959 

129.15 

277 

.2249 

96.7 

267 

.2292 

93.2 

279 

.2259 

97.4 

323 

.2  30 

112.7 

323 

.2149 

112.7 

225 

.493 

78.5 

223 

.658 

77.8 

223 

.599 

77.8 

226 

.3241 

78.9 

807 

.2592 

281.7* 

810 

.2339 

282.7* 

853 

.2773 

297.8* 

853 

.2631 

297.8* 

850 

.2917 

296.7* 

854 

.267 

298.1* 

921 

.3022 

321.49* 
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